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In the two-Higgs-doublet model (THDM), generalized- CP transformations (<fi — s> Xijtp* where 
X is unitary) and unitary Higgs-family transformations (ifi — > Uijtpj) have recently been examined 
in a series of papers. In terms of gauge-invariant bilinear functions of the Higgs fields tpi, the 
Higgs-family transformations and the generalized-CP transformations possess a simple geometric 
description. Namely, these transformations correspond in the space of scalar-field bilinears to proper 
and improper rotations, respectively. In this formalism, recent results relating generalized CP 
transformations with Higgs-family transformations have a clear geometric interpretation. 



INTRODUCTION 



The Standard Model (SM) of particle physics provides an extremely successful framework for describing the prop- 
<^> erties of the fundamental particles and their interactions. No statistically significant deviation from SM predictions 
I 'has yet been observed in collider experiments 0. Nevertheless, the experimental exploration of the scalar sector of 
^ 'particle physics is still in its infancy. The SM contains one complex doublet, hypercharge-one multiplet of scalar Higgs 
fields. But it is by no means excluded that the Higgs sector is larger than that of the SM [2]. In particular, there 
are a number of theoretical arguments suggesting a richer Higgs sector than that of the SM. For instance, the two- 
Higgs-doublet model (THDM) is attractive since it provides a viable framework for spontaneous CP violation 0, 0] ■ 
The Peccei-Quinn symmetry, originally introduced in order to solve the so-called strong CP problem [5', b], requires 
an enlarged Higgs sector and can be accommodated in the THDM. Typically, supersymmetric models require at least 
two Higgs-doublet fields [7-ll| in order to cancel potential gauge anomalies due to higgsino superpartners. Thus, 



X 



there is ample motivation to study the simplest two-Higgs-doublet extension of the SM scalar sector. 

Indeed, the theory and phenomenology of the THDM has been extensively analyzed; see 0, [lJS] and references 
therein. Among these studies of the THDM, one can find two lines of approaches. The traditional approach works 
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2 SCALAR FIELD BILINEARS IN THE THDM 



directly with the Higgs-doublet fields Ej 12-3(1. In contrast, there is a second approach that emphasizes the role of 



gauge-invariant scalar field bilinears |37H49j . A systematic use of the scalar field bilinears for the study of the stability 



and the structure of electroweak symmetry breaking (EWSB) in the THDM was first carried out in 38|, |40(. Then 



independently in [42j the scalar field bilinears were introduced and employed for the study of the CP properties of 
the THDM. This latter approach was revisited in 43 1. In a recent paper [f^], three of the present authors mainly 
employed the traditional approach in the study of Higgs-family and generalized-CP (GCP) symmetries of the THDM. 
It is the purpose of this work to present a compelling geometrical interpretation of Higgs-family and GCP symmetries 
using the formalism of scalar field bilinears. This formalism provides a powerful geometric framework that yields new 
insights and clarifies the relations among the different possible symmetry transformations of the THDM. 

In Section 2, we review the formalism of scalar field bilinears and their applications in the analysis of the THDM. 
In Section 3, we introduce the generalized-CP (GCP) transformations. These transformations are initially defined as 
transformations of scalar fields. We then revisit the GCP transformations in the formalism of scalar field bilinears 
and provide a geometric interpretation. Two useful theorems involving GCP transformations are proven in Section 4. 
Higgs-family and GCP symmetries are examined in detail in Section 5. The classification of all possible THDM 
symmetries is established in the formalism of scalar field bilinears. The constraints on the scalar potential param- 
eters due to the various possible symmetry classes is provided in a form that is covariant with respect to arbitrary 
transformations of the basis for the scalar fields. A distinction between parameter constraints derived in an arbitrary 
basis and in a specific basis is examined and clarified. Conclusions are give in Section 6. Details on the structure of 
3x3 rotation matrices that are useful for the computations performed in this paper are provided in Appendix A. 



2. SCALAR FIELD BILINEARS IN THE THDM 

The scalar sector of the THDM contains two complex doublet, hypercharge-one Higgs fields, with respect to the 
SU(2)xU(l) electroweak gauge group, denoted by 



tpi(x) 



yd*), 



with i = 1, 2 . 



(2.1) 



The tree- level THDM Lagrangian contains the kinetic term and the potential V((fi,(p2) for the Higgs fields, which 
is gauge invariant and renormalizable. These requirements imply that the scalar potential V(ipi,if2) consists only of 
quadratic and quartic terms in the fields. The conventional parametrization in the field approach reads 



V = mlxitpltpi) + m| 2 (<4v?2) - m 2 12 ((p\ip 2 ) - {ml 2 )*{<pl<pi) 
+ WiVaXvta) + W^iA^f + As(^i) 2 ] 



(2.2) 



with m\i, rn| 2 , Ai.2,3.4 real and m\ 2 , A5,6,7 potentially complex. The parameters of the scalar potential must be chosen 
such that the potential is bounded from below, is stable, and leads to the correct form for EWSB [which preserves 
U(1)em]- In addition, one may choose to impose various additional symmetry requirements on the scalar potential. 

A very convenient way to study the stability, the structure of the EWSB, and any additional symmetry requirements 
is to use scalar field bilinears 
gauge invariant bilinears as 



38 



4o[ 42 - J44I . We follow here the notation of 0, ^ and define the four independent 



Ki = ip\ip 2 



K 2 = iip^ipi - i<p\<P2, 



A", 



(2.3) 



We summarize some results from [38|, |40 1 . 
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2 SCALAR FIELD BILINEARS IN THE THDM 



• We have 

K Q > , Kl - K\ - K\ - K\ > . (2.4) 

That is, the four vectors K = (K , K) T = (K , Ki,K 2 , K 3 ) T span the forward light cone in K space. These K 
parameterize the gauge orbits of the Higgs-doublet fields. 

• A change of basis of the Higgs fields, called a Higgs- family transformation 

H - (A = u M , (2.5) 

with U = (Uij) £ U(2), corresponds to an SO(3) rotation in K space 

K -> K' = R(U)K . 

Here, the 3x3 real orthogonal matrix R(U) is obtained from 

tfa a U = R ab {U)a b . (2.7) 

Since U is continuously connected to the identity, it follows that det R{U) = 1. It is straightforward to prove 
that 

R ab {U)R bc (V) = Rac(UV) , (2.8) 

starting from eq. (|2.7[) and using the fact that the a a span the space of traceless 2x2 Hermitian matrices. Thus, 
the mapping {U , — U} i — > R(U) provides the well-known double cover of SO(3) by SU(2). An explicit formula 
for R a t,(U) is easily obtained: 

Rab(U) = ±Tr([/V[/a b ) . (2.9) 

The most general SO(3) matrix can be uniquely specified by an axis of rotation n and an angle of rotation 6 
that lies in the interval < 9 < ir. We denote the corresponding 3x3 matrix by R(n,8). The properties of 
i?(n, 9) are reviewed in Appendix A. 

• Every SO (3) rotation in K space, given by (|2.6[) with R £ SO (3), corresponds to a Higgs- family transforma- 
tion (|2.5p that is unique up to gauge transformations. 

• The most general gauge invariant and renormalizable potential V of the THDM can be written as 

V = K T £ + K T EK . (2.10) 
Here, £ and E contain the parameters which are all real, 

*-($)• ^(r £)■ < 2 -"> 

with E — E T a 3 x 3 matrix. Expressing £ and E in terms of the parameters of (|2.2p . we find 



^-2Re(m? 2 ) 

2Im(m 2 2 ) I , 77 = | 



/ Rc(A 6 + A 7 ) 
- Im(A 6 + At) 
\m\ x -m\^ J \ |(Ai-A 2 ) 



Co = |( m ll + m 22) i 

\m\ x 

'A 4 + Re(A 5 ) -Im(A 5 ) Re(A 6 - A 7 ) 

?7oo = |(Ai + A a ) + 5A3 , E = \\ -Im(A 5 ) A 4 - Re(A 5 ) -Im(A 6 -A 7 ) 

v Re(A 6 -A 7 ) -Im(A 6 -A 7 ) |(Ai + A 2 ) - A 3 , 



(2.12) 
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3 GENERALIZED-CP TRANSFORMATIONS 



• A transformation (|2.5[) corresponds to a Higgs family symmetry if and only if 

R(U)£ = £ , R(U)t] = V , R{U)ER T (U) = E . (2.13) 

In 43J basis changes of the Higgs fields as in eq. (|2.5p were considered, but with the unitary transformation U 
replaced by more general SL(2, <D) transformations. In K space this corresponds to general Lorentz transformations, 
which includes both rotations and boosts. However, the latter change the form of the kinetic terms of the THDM 
Lagrangian. Without loss of generality, one may assume that the kinetic terms of the tree-level THDM Lagrangian 
are of canonical form. Under this assumption, only unitary Higgs family transformations are permitted. 

3. GENERALIZED-CP TRANSFORMATIONS 

In this section we study the generalized-CP (GCP) transformations, 

tpi{x) ->X l3 <p*(x'), with i,je {1,2}, X = (Xy) e U(2), x=l X ),x'=l*). (3.1) 



Such GCP transformations of scalar fields have been previously considered in 5114551] 1 , and GCP transformations of 



fermions fields have also been examined in [58[. A systematic study of GCP transformations of the scalar fields of the 



THDM was carried out in [45| and [50|. The matrix X that appears in eq. (13.11) is basis-dependent. Under a change 



of basis specified by eq. (|2.5j) . the GCP transformation of eq. (|3.ip is modified to: 

xj^V) , (3-2) 

where X' is a unitary matrix given by: 

X' = UXU T . (3.3) 

As shown in [5(3], three classes of GCP transformations exist depending on the value of (GCP) 2 . Consider first 
the case of (GCP) 2 = 1, which is denoted by CP1 (sometimes called the "standard" CP symmetry transformation). 
Then, 

ipi -» Xijip* -> XijXj k ip k = ipi , (3.4) 

which implies that XX* — 1 2 (where 1 2 is the 2x2 identity matrix). Since X is unitary, the latter implies that X is 
also symmetric. Thus, eq. (|3.1[) corresponds to a CP1 transformation if and only if A is a symmetric unitary matrix. 
One can now employ the well known result that any symmetric unitary matrix X can be written as the product of a 
unitary matrix and its transpose (see e.g. Appendix D.3 of [H^] for a proof of this result). That is, one can always 
find a unitary matrix U such that X = Wll*. Performing the basis transformation given by eq. (|3.3|l then yields that 

X' = UXU T = UU\UU^)* =t 2 . (3.5) 

That is, in the case of CP1, there is always a basis choice for which X' = 1 2 . 
Next, consider the case of (GCP) 2 = —1, which is denoted by CP2. Then, 

ipi -> X^^* -> XijX* k (f k = -ipi , (3.6) 



1 For the discussion of CP transformations in the SM see for example |55l45 
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3 GENERALIZED-CP TRANSFORMATIONS 



which implies that XX* = —\^- Since X is unitary, the latter implies that X is also antisymmetric. Thus, eq. (|3.1[) 
corresponds to a CP2 transformation if and only if X is an antisymmetric unitary matrix. The most general anti- 
symmetric unitary 2x2 matrix X is 2 

X = e** € , where e = ia 2 = -e T = -er 1 = ^ ° ^ , (3.7) 

where tp is an arbitrary phase. It then follows that a unitary matrix exists, U = e _l ^/ 2 e, such that 

X' = UXU T = e. (3.8) 

That is, in the case of CP2, there is always a basis choice for which X' = e. 

Finally, we consider the case of (GCP) 2 = XX* ^ ±la, which is denoted by CP3. In this case, it is always possible 
to perform a basis transformation such that in the new basis, X is transformed into: 3 

X> = UXU^( COS0 ff Sm f), (3.9) 
\ — sin cos J 

where < 6 < ir/2. 

To prove this result, we first note that since X is a unitary matrix, det X = e 2lx is a pure phase. Following eq. (13.31) . 
we shall perform a basis transformation such that 

det X' = det(UXU T ) = 1 . (3.10) 

This can always be done provided that 

U = e- l ^ 2 U, (3.11) 
where U is an SU(2) matrix. It is convenient to define: 

X = e- ix X, (3.12) 

in which case det X = 1 and 

X' = UXU T . (3.13) 

A general SU(2) matrix U satisfies: 

U = eU*e-\ (3.14) 

where e is defined in eq. (|3.7p . Eq. p. 141) expresses the well known equivalence of the irreducible two-dimensional 
representation of SU(2) and its complex conjugate. Inserting the transpose of eq. (|3.14[) into eq. (|3.13[) yields: 

X'e = UXe C/ f . (3.15) 

It is convenient to define: 



^f 1 

V2 \i 



C/^^l 1 ! 7, (3.16) 



— % 



2 Conversely, it is straightforward to show that if X^a a X = —a a *, where X = exp(i9fi • <t/2), then X is proportional to a 2 , i.e. X is an 
antisymmetric unitary matrix. 

3 Eq. I|3.9|l is an example of a canonical form for unitary congruence. For a comprehensive mathematical treatment, see |6dl | (note in 
particular Corollary 8.7). In the physics literature, eq. 113.91 1 first appeared in 53] and was further generalized in [6lll . 
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where V is a unitary matrix (such that det V = i). Since Xe is an SU(2) matrix, it follows that the two eigenvalues of 
Xe are complex conjugates of each other, denoted below by e ±l< ^, where the real angle <f> is denned modulo n. Then, 
we can choose V to be the unitary matrix that diagonalizes Xe, 

VXeVt= ^ °]. (3.17) 



e 



Inserting eqs. (|3.16|) and (|3.17|l into eq. (|3.15|) yields 



X'= Sln t " COS H. (3.18) 



cos m sin i 



Finally, we define <j) = 6 + tt/2 to obtain the desired form given by eq. (|3.9I) . 

The angle cj> (and hence the angle 9) is defined modulo n. Thus, it is convenient to establish a convention in 
which \9\ < tt/2. However, we are free to redefine U — > cr 1 C/ ', which has the effect of changing the overall sign of 6*. 4 
Consequently, it is always possible to find a basis transformation such that X 1 takes the form given by eq. (|3.9[) . where 
< 6 < tt/2. Moreover, 

XX* = U^X'U*U T X'*U = U^ ( C0SW S ' mW )u, (3.19) 

^ sin 26» cos26»y V ' 

so that 9 = corresponds to the case of CP1 [eq. (|3.5|) ]. 9 = n/2 corresponds to the case of CP2 [eq. (|3.8j) ]. and 
< 9 < tt/2 corresponds to the case of CP3 [eq. (|3.9[) ]. 

Summarizing the above results, it follows that in a suitable basis for the scalar fields the matrix X in p.ip can 
always be brought to the form 



with < 9 < tt/2. (3.20) 



/ cos 9 sin 9 \ 
\ — sin 9 cos 9 I 

The classification of GCP symmetries established above is [5(3]: 

• CP1 if 9 = 0, 

• CP2 if 9 = n/2 and 

• CP3 if < 9 < n/2. 

We now demonstrate how the classification of GCP symmetries can be understood in the formalism of field bilinears 
employed in [45]. To make the present paper self-contained, we shall repeat some of the derivations of [45J in the 
analysis that follows. 

In the notation of eq. (|2.7[) . we define the SO(3) matrix R a t,(X) via: 

X^a a X = R ab {X)a\ (3.21) 

where X is the unitary matrix that specifies the GCP transformation [cf. eq. (|3.ip j. It is convenient to introduce the 
improper rotation matrix, 

i? 2 = diag(l,-l,l). (3.22) 



4 For any 2x2 matrix A, the matrix o^Acr 1 is related to A by an interchange of the two diagonal elements and an interchange of the two 
off-diagonal elements. 
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The matrix R 2 describes the reflection through the 1-3 plane in K space. Similarly, we introduce R\ and R3 as the 
reflections through the 2-3 and 1-2 planes, respectively, 5 

Ri = diag(-l,l,l), 

_ bv > ) n (3 23) 
fl 3 = diag(l,l,-l). 

In particular, note that 

a a * = a aT = (R 2 ) ab a b . (3.24) 

The scalar field bilinears of eq. (|2.3p can be rewritten as: 

K, = <p\o» j(Pi , (3.25) 

where cr M = (1,ct). Then, the GCP transformation defined in eq. (|3.ip corresponds to 

*m - <4<-^ -> (X^lY^X^ = pJ W (^) M . (3.26) 

If /x = 0, then Xt(T J£ = X f X = 1 2 , and so K -> AT (where we suppress the coordinates x and x'). If /i = a = 1, 2, 3, 
then one may use eqs. 13.21(1 and (|3.24|) to obtain: 

K a ^<p\<p k (X^a a X) M 
= (pl(p k Rab(X)a b ke 
= ip\w k R ab {X)(R-i)bcGtk 

= RacK c , (3.27) 



where R is the improper rotation matrix: 



That is, 



R=R(X)R 2 . (3.28) 



K (x) -> A' (a;'), 
K(i)->SK(i') 



(3.29) 



which reproduces the result obtained in section 3 of 45 1. 

Under each of the three classes of GCP transformations, the improper rotation matrix R = R(X)R 2 must satisfy an 
appropriate constraint equation. To derive the relevant constraint, we start with the complex conjugate of eq. (|3.2ip . 
Using eq. (|3.24[) . it then follows that 

X T (R 2 ) ab a b X* - R ab (X)(R 2 ) bc a c . (3.30) 

Employing eq. (|3.21[) once more yields 

(R 2 )abRbc(X*)a c = R ab (X)(R 2 ) bc a c . (3.31) 
Since the a a are linearly independent and span the space of traceless 2x2 Hermitian matrices, eq. (|3.31[) yields: 

R(X*) = R 2 R(X)R 2 , (3.32) 



5 Here and in the following R, R a , . . . and _R, R a ,... denote proper and improper rotation matrices with determinant +1 and — 1, 
respectively. 
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3 GENERALIZED-CP TRANSFORMATIONS 



after using (R2) 2 = I3 (where I3 is the 3x3 identity matrix). Finally, we multiply eq. (|3.32l) on the left by R{X) 
and make use of eq. (|2.8[) to obtain: 

R{XX*)=Tf, (3.33) 

where the improper rotation matrix R = R(X)R 2 was introduced in eq. (|3.28l) . 

Consider separately the cases CP2 (where X is antisymmetric and XX* = — 12) and CPl (where X is symmetric 
and XX* = 1 2 ). In both cases, eq. (j3~2"Tl) yields R(XX*) = R(±t 2 ) = 1 3 , and it follows that T? = 1 3 . That is, 
R is either a reflection matrix corresponding to a reflection through some plane in K space or an inversion (or point 
reflection) through the origin in K space. In eq. (|3.7[) . we noted the most general form for X in the case of CP2 is 
given by X — e^ia 2 , where ip is an arbitrary phase. Inserting this result into eq. (|3.21[) and making use of eq. (|3.24j) 
yields: 

R ab (X)a b = a 2 a a a 2 = -a a * = -{R 2 ) ab a b , (3.34) 

and we conclude that R(X) = — R%, Hence, eq. (|3.28|) yields R — — 13, which corresponds to an inversion (i.e., a point 
reflection through the origin in K space). This case is a GP^ transformation in the notation of 45]. In contrast, 



2 

for the case of CPl, where X is a symmetric unitary matrix, we have R = I3 and R 7^ — 13. An improper rotation 
matrix of this type corresponds to reflections through planes in K space. In particular, an SO(3) matrix R exists 
such that 

R = RR 2 R T . (3.35) 

To prove eq. (|3.35l) . simply choose a basis where X = 1 2 in which case R(X) — 1 3 and R — R 2 . Then, rotate in 

~ 2 

K space to an arbitrary basis using the rotation matrix R to obtain eq. (|3 .35|) . One can easily check that R = I3 



and R 7^ — 13 as required. This case corresponds to a CPg"' ) transformation in the notation of [45j. In summary, for 
those GCP transformations whose square is equal to the unit transformation when acting on t he g auge invariant field 
bilinears, we must have R 2 = I3. The resulting classification of jif^ is then related to that of [50] as follows: 

• CP^"': reflections through planes in K space, R — RR 2 R T with R e SO(3) CPl transformations. 

• CP^: a point reflection through the origin in K space, R——I3 <^=> CP2 transformations. 

The case of CP3 transformations was not considered in detail in 45]. In this case, XX* 7^ ±1.2, which implies that 

2 

R(XX*) 7^ I3. Hence, eq. (|3.33|) yields R 7^ I3, which means that the improper rotation R is not a reflection matrix 
or an inversion. To obtain an explicit form for R = R(X)R 2 , it is convenient to choose a suitable basis in which X is 
given by eq. (|3.20p . which can be rewritten as: 

X = l 2 cos6> + icr 2 sm<9. (3.36) 
We can use eq. (|2.9I) to obtain the corresponding rotation matrix R a i,(X). Evaluating the relevant traces, one obtains: 

Rab(X) = S ab cos 26 + 25 a2 5 b2 sin 2 6 + e ab2 sin 29 . (3.37) 
Hence, in a basis in which X is given by eq. (|3.20[) . 

/cos 29 -sin26>\ 
R=R(X)R 2 =0 -1 (3.38) 

^sin26> cos26»y 

As expected, for a CPl transformation 9 = and R = R 2 [in a suitable basis in which X = 1 2 ], and for a CP2 

transformation 9 = ir/2 and R = — 13. The case of < 9 < tt/2 corresponds to a CP3 transformation, in which R is 

2 

an improper rotation matrix that is not a simple reflection or inversion (i.e. R 7^ I3). As in eq. (|3.35|) . the form for 
R in a general basis is related to eq. (|3.38p by an orthogonal similarity transformation, 

R = RR(X)R 2 R T , (3.39) 
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4 TWO THEOREMS INVOLVING GCP TRANSFORMATIONS 



for some SO (3) matrix R, where R(X) is given by eq. (13.37[) . 

Thus, we have reproduced above the result proved in [451 ] . Every transformation of the field bilinears given 
by eq. (3.29)) . where R is any improper rotation matrix, corresponds to a GCP transformation of the fields as specified 
in eq. (3.1)) . This field transformation is uniquely determined by R up to gauge transformations. In analogy with 
eq. (|2.13p , a GCP transformation of the form given by eq. (I3.29[) corresponds to a GCP symmetry if and only if 

R£ = £ , RV = V, RElf = E . (3.40) 



TWO THEOREMS INVOLVING GCP TRANSFORMATIONS 



In [50| , it was suggested that all symmetries of the THDM could be expressed in terms of products of GCP symme- 
tries. In this section, we prove two simple theorems that demonstrate that all Higgs family and GCP transformations 
can be expressed in terms of products of CP1 transformations. 

Theorem 1: Any Higgs-family transformation can be considered as a product of two CP1 transformations. 

Theorem 2: Any GCP transformation is either a CP1 transformation or a product of three CP1 transformations. 

These theorems are new. Some particular cases were considered in [501 ] but even there they referred only to symmetries 
where the form of the potential could play a role. The results presented here apply to the more fundamental trans- 
formations themselves. 

Proof of Theorem 1. Let X and Y be two arbitrary symmetric unitary matrices. Consider the product of the 
corresponding CP1 transformations, 

tp^x) = X jk tp* k (x') , (4.1) 
tf{x) = Y l3V '*{x') = Y l3 X* kVk {x) = U ik tp H {x) , (4.2) 

where U = Y X* is a unitary matrix. The theorem is proven if we can show that an arbitrary unitary matrix U is the 
product of two symmetric unitary matrices. But this last statement is easy to prove. First, we diagonalize U with a 
unitary matrix W, 

U = WDW^ , (4.3) 

where D is a diagonal matrix of phases corresponding to the eigenvalues of U. Then define the following two symmetric 
unitary matrices: 

Si = WW T , S 2 = W*DW J{ . (4.4) 



It immediately follows that 



= WW T W*DW 1 = WDW^ = U , (4.5) 



which shows that any unitary matrix can be written as a product of two symmetric unitary matrices. 

It is also instructive to prove Theorem 1 in K space. Consider an arbitrary Higgs-family transformation [cf. eqs. (|2.5[) 
and (|2.6[) ]. Every proper rotation matrix R is a rotation about an axis and can be represented, in a suitable basis, as 



R a = 



cos a — sin a 
sin a cos a 
:i 



with < a < 7r. (4.6) 
This rotation can also be generated by two reflections, R a /2 and i? Q , as illustrated in Fig. [TJ 
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second reflection plane first reflection plane 

\ \ K 3 / 

\ / 

\ K / 

\ f 7 

\ / \/ 




FIG. 1: The reflections R a /2 and R a illustrated in the K1-K2 plane. The reflection R a /2 brings the arbitrary vector K to K. 
Then, R a brings K to K'. 

The first reflection, R a /2, is through the plane containing the K3 axis and the line of the angle a/2 in the K\-K 2 
plane. The second reflection, R a , is through the plane containing the K3 axis and the line of the angle a itself. 
The arbitrary vector K goes into K under the reflection R a /2 and further into K' under R a . In Fig. Q]we show 
the projections of these vectors onto the K\-K 2 plane. Since the third component of the vectors are not changed 
under R a /2 and R a we see that K' gives exactly the vector rotated by the angle a, that is, K' = R a K, where R a is 
given by eq. (|4.6p . To see all this formally we consider Rp as in eq. (|4.6p but with a replaced by an arbitrary angle /3 
and define 

sin 2/3 

Rp = RpR 2 Rj = I sin 2/3 -cos 2/3 0|. (4.7) 

Note that (Rp) 2 = I3, which indicates that Rp is a pure reflection. In particular, Rp describes the reflection through 
the plane containing the K 3 axis and the line of angle /3 in the K\-K 2 plane. We have Rj — R-p, R 2 RpR 2 = R~p, 
and Rp 1 Rp 2 = Rp 1 +p 2 . It is then easy to see that we get R a from eq. (|4.6[) via two reflections, 

This is the K space equivalence of the statement that any Higgs family transformation is equivalent to the product 
of two CP1 transformations. 

Proof of Theorem 2. If we perform three successive CP1 transformations of the form given by eq. (|3.1|> . with X 
given by Xi , X 2 and X 3 , respectively, then the resulting transformation is of the form of eq. (|3.ip with 

X = X 3 X;Xi , (4.9) 




where the Xi are symmetric unitary matrices. However, any unitary matrix X can be written in the form of eq. (|4.9[) . 
This follows from the fact that any unitary matrix can be written as the product of two symmetric unitary matrices, as 
proven in eq. (|4.5j) . Thus, one can simply choose one of the matrices Xi in eq. (|4.9j) to be the identity matrix 6 (which 



6 Of course, this does not imply that a GCP transformation can be expressed as the product of two CP1 transformations. Each CP1 
transformation involves the complex conjugation of the scalar fields, so one requires a product of odd number of CP1 transformations to 
express a GCP transformation. The product of two CP1 transformations is a Higgs family transformation as demonstrated in Theorem 1. 
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is of course a symmetric unitary matrix). Finally, we note that the product of an odd number of GCP transformations 
is a GCP transformation, as this follows trivially from eq. (|3.1[) . The proof of Theorem 2 is then complete. 

Again, it is instructive to prove Theorem 2 in K space. Consider an arbitrary GCP transformation specified by 
eq. (|3.29[) . In a suitable basis R has the form in eq. (|3.38[) . For 9 — we are finished, since R = R2 which corresponds 
to a CP1 transformation. For 9 ^ we make a basis transformation in eq. (|3.38[) exchanging the 2 and the 3 axes, 
which yields an improper rotation, 



R 



20 



( cos 26* - sin 29 

sin 29 cos 2(9 | . (4.10) 
y 



We can represent this matrix as a product of three CP1 transformations by multiplying eq. (|4.8p by -R3 and replacing 
a by 29. Thus we obtain 

RzRieRe — (i?i?2^ T ) (-R20-R2 ^29) (R9R2RJ) = R26 > (4-11) 

where 



R=\0 -1 . (4.12) 




Eq. (|4.11[) is the K space equivalence of the statement that an arbitrary GCP transformation either is a CP1 trans- 
formation or the product of three CP1 transformations. 

Note that Theorem 2 is equivalent to the statement that an arbitrary improper 3x3 rotation matrix can be 
expressed as a product of three reflection matrices. This statement can be proved directly. First, we note that any 
proper 3x3 rotation matrix can be written as the product of two reflection matrices as in eq. (|4.8[) . However, 
any improper rotation matrix can be written as the product of a reflection matrix and a proper rotation matrix. 
Combining these two statements yields the desired result. 7 

5. HIGGS-FAMILY AND GCP SYMMETRIES 

In Sections 3 and 4, we analyzed Higgs-family and GCP transformations. In this section, we shall address the 
corresponding symmetries. In the scalar sector of the THDM, there are six nontrivial inequivalent classes of symmetries 
that can be exhibited by the tree-level THDM scalar potential. The complete classification of the THDM symmetry 
classes, which are listed in Table I, was first given in 44] and subsequently analyzed in [50I ] . 

The following points are noteworthy. First, consider the case in which a Higgs-family or a GCP transformation is a 
symmetry transformation of the THDM. The representation of such a symmetry transformation as a product of other 
transformations, (e.g. a product of CP1 transformations as discussed in Section 4), does not automatically imply that 
the individual factors of the product are also symmetries of the THDM. Second, the analysis of this paper is confined 
to the scalar sector of the THDM. Ultimately, one must also include the Yukawa couplings to fermions in the theory 
in the discussion of symmetry transformations. The structure of the symmetries of the scalar sector may not be 
respected by the Yukawa sector. As an example, consider a field transformation Si that can be written as a product 
of two others, S\ = S2S3. Suppose that by imposing the symmetry Si, the symmetries S2 and S3 are automatically 
respected by the scalar sector. This property is not guaranteed to hold for the Yukawa sector. In particular, it is 
possible that imposing S2 and S3 separately as symmetries may lead to stronger restrictions as compared with the 



imposition of the symmetry Si by itself. This is indeed the case in the model which was studied in 4y, |47|, |49| . 



7 In particular, the inversion matrix — 13 can be written as the product of three reflection matrices. Starting from eq. I|4.8j l. we note that 
—Ra is a reflection matrix, in which case, —RaRaR a /2 = — 13- 



11 



5.1 Symmetries and single rotations in bilinear space 



5 HIG GS-FAMILY AND GCP SYMMETRIES 



symmetry class 


rotation matrices in a generic basis 


rotation matrices as a product of CP1 reflection matrices 


Z 2 


R(h, n) 


R-k = R2R1, for n = z [see eq. (15.41)] 


U(l) 


R{n,0), (0<6»<tt) 


R?.a = Ro.aRa. for n = z [see ea. (|5.23|l] 


SO(3) 


±R(hi,9) , ±R(n 2 ,G) , (ni x n 2 + 0) 


See section [5731 


CP1 


-R(n,n) 


Ro. for h = v [see ea. (|5.34l)] 


CP2 


-R(n,Q) = -1 3 


-1 3 = R3R2R1 [see eq. (|5.44|)] 


CP3 


-R{h,6), (0<6»<tt) 


H 2 0i = R3R2e 1 Re 1 , for n = z [see eq. (|5.46|l] 



TABLE I: The symmetry classes and the corresponding proper [and improper] rotation matrices R [and R] that generate 
the symmetry classes in K space. A general rotation matrix R(n, 9) is uniquely determined by an axis of rotation n and 
rotation angle 8 (see Appendix A). The identity class generated by R(n,0) = I3 is trivial and is not explicitly displayed. The 
decomposition of R and R in terms of products of CP1 reflection matrices is also given in a particular basis. In the case of the 
CP2 and SO(3) symmetry classes, the corresponding rotation matrices are invariant with respect to basis transformations. The 
generation of the SO (3) symmetry class requires two rotation matrices (each of which may be a proper or improper rotation 
depending on the overall sign choice), where fi\ and n 2 are non-collinear. 



5.1. Symmetries and single rotations in bilinear space 

We shall now use Theorems 1 and 2 of Section 4 to obtain a new derivation of the statement (given in [5^]) that 
all possible symmetries of the scalar sector of the THDM can be reduced to multiple applications of the standard CP 
symmetry in suitable bases. In Table Q] we list the possible classes of THDM symmetries, along with the equivalent 
result in K space, where the Higgs family and generalized CP transformations are generated by proper or improper 
rotations, R or R, respectively. Also shown is the decomposition of R and R in terms of products of CP1 reflection 
matrices in a particular basis choice. 8 Note that it is sufficient to require the invariance under a single but suitable 
rotation in order to generate each symmetry class, with the exception of the SO(3) case. 9 

The transformation matrices R and R generating the respective symmetry classes can be written as products of 
CP1 transformations using Theorems 1 and 2 of Section 4. Of course, one still needs to check if the imposition of the 
CP1 factors as symmetries is equivalent or more restrictive than the corresponding proper or improper rotation, R 
or R, alone. Let us now derive and discuss the results shown in Table U in detail. 

• Z2 symmetry 

The Z2 symmetry 62, 6^ . 



(5.1) 



<pi(x) -t <pi(x), 

ipi{x) -> -<f2{x), 

corresponds in K space to a rotation by 7r around the third axis, 

K(x) -> R^Kix) , (5.2) 

where R^ is as in eq. (|4. 6[) with a — n. That is, 

'-1 0^ 

/?■_, = lh = I 0-1 I . (5.3) 

1, 



8 In the CP2 symmetry class, the rotation matrices are independent of n, so that the basis choice in this case is fixed by the direction of 
an eigenvector of E corresponding to one of its non-degenerate eigenvalues. 

9 For the SO(3) symmetry class, if a basis for the scalar fields is appropriately chosen, then it is possible to generate the SO(3) symmetry 
class via a single rotation matrix, as shown in Section 5.3. 



12 



5.1 Symmetries and single rotations in bilinear space 



5 HIG GS-FAMILY AND GCP SYMMETRIES 



This result is most easily obtained by taking U = a 3 in eq. (|2 .9[) . 
The application of Theorem 1 gives the decomposition [cf. eq. (|4.8j) ]: 

R-rr = RttRjt/2 = R2R1 ■ (5.4) 

Requiring Z2 to be a symmetry means that the THDM parameters specified in eq. (|2 . 12[) must satisfy eq. (|2.13l) with 
R{U) replaced by R^, which yields: 

6 =6 = 0, m=V2 = 0, E 13 = E 23 = 0. (5.5) 

On the other hand, requiring the CP1 transformations corresponding to both R\ and R2 to be separate symmetries 
gives, in addition to eq. (|5.5[) . 

E 12 = 0. (5.6) 

Thus, imposing the CP1 symmetries R± and R 2 yields a stronger constraint than the Z2 symmetry alone. Nevertheless, 
the scalar potential subject to the CP1 symmetries R\ and R 2 is physically equivalent to the scalar potential subject 
to the Z2 symmetry, since the two scalar potentials are relate by a change of basis. This is easily proved by performing 
a change of basis characterized by an SO(3) matrix R [cf. eq. (|2.6I) ] with R\ 3 = R 23 = R 3 i = R 32 = and R 33 = 1. 
Note that £ and rj are invariant under this change of basis, whereas E — > RER T . One can choose R such that E is 
diagonal, which confirms that the two scalar potentials are related by a basis transformation. 
It is also instructive to introduce the permutation symmetry II2, 

<pi(x) -> ip 2 (x), 

(5.7) 

<p 2 (x) -> <pi(x). 

In fact, the Z2-symmetric scalar potential and the ^-symmetric scalar potential are related by a basis transforma- 
tion 25, 5(3] (and are hence physically equivalent). To obtain the K space description of II2, simply insert U = a 1 
into eq. ()2.9|) . which yields: 




Rn 2 = I 0-1 I . (5.8) 

Requiring II2 to be a symmetry implies that the THDM parameters specified in eq. (|2. 12|) must satisfy eq. (|2.13[) with 
R{U) replaced by i?n 2 > which yields: 

6=6 = 0, ry 2 = % = 0, E 12 = E 13 =0. (5.9) 



It is possible (although not particularly illuminating) to construct the basis transformation that relates eqs. (|5.5p 
and (|5.9[) . However, it is more useful to identify the most general transformation in K space that corresponds to the 
presence of the Z2 symmetry specified by eq. (|5.1|) in some basis. This can be accomplished by starting in a basis 
where eq. (|5.2j) is satisfied and transforming to an arbitrary basis. 

Under a basis transformation specified by the U(2) matrix U [cf. eqs. (|2.5[) and (|2.6p ]. we define R = R(U), where 
R(U) is given by eq. (|2.9[) . Then, the Z2 symmetry transformation, K(x) — > i? 7r K(x), is transformed to 

RK(x) RR{z, njipRKix) , (5.10) 

where R^ = R(z,ir) is a rotation by 180° about the z-axis, and the expression R T R = I3 has been conveniently 
inserted. Using eq. (|A.20j) given in Appendix A, it follows that 

K'(x) = RK(x) = R(n,ir)K'(x) , n = Rz, (5.11) 

which is the form of the Z2 symmetry in the new basis. 
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Henceforth, we drop the primed superscripts. The most general transformation in K space that corresponds to the 
presence of the Z2 symmetry in some basis is given by: 

K(x)^ R{n,ir)K(x), (5.12) 

where R(n,Tv) is a rotation by 180° about an axis that is parallel to the unit vector h. As noted in Appendix A, 
R(n,ir) possesses one non-degenerate eigenvalue equal to +1 and two degenerate eigenvalues —1. The eigenvector 
corresponding to the non-degenerate eigenvalue +1 is identified as the rotation axis ra, since 

R(h,Tr)h = h (5.13) 

is just the geometrical statement that the rotation axis is unaffected by the rotation. Under the symmetry governed 
by eq. (|5.12j) . the THDM parameters specified in eq. (|2.12[) must satisfy eq. (|2.13[) with R(U) replaced by R(n,Tv), 
which yields: 

£ and r] are eigenvectors of R(h,%) with eigenvalue +1 <^=^ £ and r\ are parallel to h, (5-14) 

ER(h, tt) = R{h, ir)E Eh is parallel to h . (5.15) 

Eq. (|5.14p is a consequence of the fact the eigenvalue +1 of R(h, tt) is non-degenerate, which implies that any vector 
v that satisfies R(h,TT)v = v must be proportional to h. To derive eq. (I5.15|) . we note that h is a simultaneous 
eigenvector of R(fi,w) and E. In particular, R(h,n)Eh — ER(h,Tr)h = Eh, where the last step follows from 
eq. (|5.13|) . Hence either Eh = or Eh is an eigenvector of R(h, tt) with eigenvalue +1. Since the latter is non- 
degenerate, it follows that 10 

EijHj oc ni. (5.16) 

One can easily check that eqs. (|5.14[) and (|5.16l) reduce to eq. (|5.5p or eq. (|5.9p when h = z or x, respectively. 

Eq. (|5.15|) implies that the eigenvectors of E can be chosen to be simultaneous eigenvectors of R(h,ir). Since E 
is a real symmetric matrix, these eigenvectors can be chosen to be orthonormal. We denote these eigenvectors by 
{h , rh , h X rh}. We have already noted that h is an eigenvector of E by virtue of eq. (|5 . 1 6[) . Thus, the other two 
eigenvectors of E must satisfy: 11 

R(h,ir)rh = —rh , R(h, 7r)(n X rh) — —h X rh . (5-17) 

Because of the two-fold degeneracy of the eigenvalue —1 of R(h, tt), it is possible to perform orthogonal transformations 
within the two-dimensional subspace spanned by rh and h X rh that leave the form of the Z2 symmetry transformation 
given in eq. (|5.12p unchanged. This simply means that the form of the 1i symmetry transformation does not uniquely 
specify the basis in K space. To fix the basis completely, one must specify h and the eigenvectors of R(h,n) 
corresponding to the two- fold degenerate eigenvalue — 1. 

In summary, the Z2 symmetry corresponds to K(x) — > i?(n, 7r)K(x) for some choice of h. Imposing this symmetry 
on the scalar potential requires that 

£ and rj are parallel to h , (5.18) 
E is diagonal with respect to the basis {h , rh , h X rh} , (5.19) 

where rh is a simultaneous eigenvector of R(h, tt) and E, with R(h, Tr)rh = —rh. In this case, the choice of h and rh 
uniquely fixes the basis in K space. 



10 One can derive eq. 15.1611 more explicitly by employing eq. J A.6ft in ER(n,Tr) = R{n,ir)E . It follows that niUjEj^ = n^njEji. Using 
the fact that E T = E, one obtains EijUj = Crii, where the constant of proportionality is identified as C = rijEj^n^. 

11 In general, the eigenvectors rh\ and rn.2 defined in eq. l lA.151 1 are not expected to be eigenvectors of E. In this case, rh is some linear 
combination of mi and TT12, and similarly for n X rh. In particular, E is not generally diagonal with respect to the basis {h , rh\ , rh.2}. 
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• U(l) Peccei-Quinn symmetry 

The U(l) Peccei-Quinn symmetry 0,0,0] requires invariance under 

<pi(x) -> e~ ie ipi(x), 
ip 2 (x) -> e l8 (p 2 (x) , 



(5.20) 



with an arbitrary angle 8, which is defined modulo tt. By taking U = \ 2 cos 9 — ia 3 sind in eq. (|2.9p . we obtain the 
corresponding symmetry transformation in K space, 

K(x) -> R 2e K(x) , (5.21) 

where R 2 g is as in eq. (|4.6p with a = 28, with < 8 < tt. That is, 

/cos 28 -sin 28 
sin 28 cos 26» | . (5.22) 

V 



-Ru(l) — 




The application of Theorem 1 gives the decomposition [cf. eq. (14.8[) 

R 2 g = R 2 gRg. (5.23) 

Requiring the U(l) Peccei-Quinn transformation to be a symmetry implies that the THDM parameters specified in 
eq. (|2.12|) must satisfy eq. p.!3|) with R(U) replaced by R 2 g (for all possible values of 8), which yields: 

6=6 = 0, 771=772 = 0, J5 = diag(/ii,/ii,/i 3 ). (5.24) 

It is straightforward to check that imposing the two CP1 symmetries Rg and R 2 g separately is equivalent to requiring 
invariance under R 2 g. 

Remarkably, it is sufficient to require invariance of the scalar potential given by eq. (|2.2j) under R 2 g for any single 
particular value of 28q ^ (mod tt) 64]. 12 That is, invariance under R 2 g for any 28q ^ (mod tt) implies eq. (|5.24|) . 
which in turn implies invariance under R 2 g for all values of 8. 

More generally, we consider the possibility that a basis transformation is required to identify the symmetry specified 
in eq. (I5.20[) . First, we shall rename 28 by 8 and rewrite eq. (|5.21j) as 

K(x) -> R(±z,8)K(x), (5.25) 

where < 8 < tt. Note that we must allow for both signs of ±z in order to cover the entire U(l) Peccei-Quinn 
group manifold. Under a basis transformation specified by the U(2) matrix U [cf. eqs. (|2.5j) and (|2.6j) ]. we define 
R = R(U), where R(U) is given by eq. (12.91) . Then, with assistance from eq. (|A.20[) . the U(l) Peccei-Quinn symmetry 
transformation, eq. ()5.25|) . is transformed to 

K'(x) -> RR(±z, 8)R T K\x) = R(±ii, 8)K'(x) , (5.26) 

where n = Rz and i?(n, 8) is a rotation by 8 about the axis n. 

Thus, dropping the primed subscripts, the most general transformation in K space that corresponds to the presence 
of the U(l) Peccei-Quinn symmetry in some basis is given by: 

K(x) -> R(n,6)K(x) , for0<6><7r. (5.27) 



Note that if 9q = it/2 then R2e = Rtt, which generates the Z2 symmetry class treated previously. What is perhaps more surprising is 
that if Go = Tr/n for any integer n > 2, then the invariance of the scalar potential under i?2e [which generates a Z n subgroup of U(l)] 
implies invariance under the full U(l) group. This latter result is a consequence of the fact that the scalar potential of eq. 112.21 1 contains 
no terms of dimension greater than four. If one relaxes this condition, then new symmetry classes can arise that are associated with 
discrete subgroups of U(l) of order d > 2. 
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Note that we have excluded the case of 9 = 0, which corresponds to the identity transformation, and the case of 9 = 7T, 
which corresponds to the Z2 symmetry transformation treated previously. When 6> 7^ (mod ir), R(h,9) possesses 
three non-degenerate eigenvalues: +1, e ie and e~ lB , and ft is the normalized eigenvector of R(h,9) with eigenvalue 
+1. It is convenient to introduce normalized eigenvectors rh and rh*, corresponding to the eigenvalues e and e~ te , 
respectively. For further details, see Appendix A. 

Under the symmetry governed by eq. ()5.27[) . the THDM parameters specified in eq. (|2.12[) must satisfy eq. (|2.13j) 
with R(U) replaced by R(n,9), which yields: 

£ and r] are eigenvectors of R(n, 9) with eigenvalue +1 •<=>■ £ and rj are parallel to n, (5.28) 

ER(n, 9) — R(n, 9)E Ev is parallel to v for v = h , rh , rh* . (5.29) 

The derivation of eq. (|5.29|) is similar to the one given in the case of the 1 2 symmetry above. We thereby obtain 
two conditions on the matrix E. The first condition, 

EijUj cx rii , (5.30) 

coincides with eq. (|5.16[) . To derive the second condition, we note that R(h,9)Erh = ER(h,9)rh = e Erh, which 
implies that either Erh = or Erh is an eigenvector of R(n,9) with eigenvalue e l6 '. Since 9 7^ (mod n) by 
assumption, the latter is non-degenerate, and it follows that 13 

Eijrrij oc m-i. (5.31) 

This equation and its complex conjugate imply that rh and rh* are eigenvectors of E, whose eigenvalues are complex 
conjugates of each other. But E is a real symmetric matrix, which implies that its eigenvalues are real. It follows 
that there is (at least) a two-fold degeneracy among the eigenvalues of E. 

As noted below eq. (|A.14[) . the eigenvector rh is independent of the rotation angle 6, assuming that 9 7^ (mod ir). 
Consequently, the constraints on the scalar potential parameters [governed by eqs. (|5.28|) , (|5.30p and (|5.31|) ] do not 
depend on 9. That is, the invariance of the scalar potential under R(n, 9) for any single particular value of 9 7^ 
(mod 7r) yields the U(l) Peccei-Quinn symmetry, which in turn implies the invariance of the scalar potential under 
R(h, 9) for all values of 9. 

One can check that eqs. (|5.28[) . (|5.30|) and (|5.31|) reduce to eq. (I5.24p in the basis where h = z and rh = -^={x — iy). 
For example, after inserting rh into eq. (|5.31[) . and taking into account that E is a real symmetric matrix, we obtain: 

E n - iEu = %Ei3 + E 22 , E 13 - iE 23 = . (5.32) 

Taking the real and imaginary parts of the above equations, it follows that E is a diagonal matrix with En = E 22 as 
stated in eq. (|5.24p . Indeed, the above computation is valid for any choice of 9 7^ (mod 7r), as noted above. Due 
to eqs. (|5.30p and (|5.3ip . the vectors comprising the orthonormal set {n,rh,rh*} are simultaneous eigenvectors of 
R(n, 9) and E. This means that E is diagonal with respect to the basis {h , rh , rh*}. Since rh and rh* are uniquely 
determined by h [cf. eqs. (|A.13|) and (|A.14[) ] when 9 7^ (mod ir), it follows that h uniquely fixes the basis in K space. 

• CP1 symmetry 

In Section 3, we showed that one can always find a basis in which the CP1 transformation in K space is given by 

K(a;) -> R 2 K(x') , (5.33) 



One can derive eq. I|5.31|l more explicitly by employing eq. l)A.5| l in ER(n, 9) = R(h, 9)E. This yields two constrains: riirijEjk = n^njEji 
(also obtained in the case of the Z2 symmetry, as noted in footnote llUt and e i j^_E jfe n^ = ej^iEijUf. If we multiply the latter equation by 
rrik and employ E T = E, En oc n and n-rh = 0, it then follows that Eijrrij = rem,, where the constant of proportionality is identified 
as k = ^(Tr E — rii Ey nj ) . 
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That is, 

- - I' ° °\ 

Rcpi = Ri = -1 . (5.34) 
V° V 

Requiring CP1 to be a symmetry implies that the THDM parameters specified in eq. (12.12)) must satisfy eq. p.40j) 
with R = i?cpi, which yields: 

6 = , 772 = , E 12 = E 23 = 0. (5.35) 

As expected, eq. (|5.35[) is equivalent to the statement that in a basis in which the CP1 transformation is given by 
eq. (|5.33p , all the parameters of the scalar potential specified in eq. (|2.2[) are real. 

As in the case of the Z2 symmetry, there still remains some freedom to perform a basis transformation while 
maintaining the form of the CP1 symmetry given in eq. (|5.33p . By performing a change of basis characterized by an 
SO (3) matrix R with R12 = R23 = R21 = R32 = and R22 = 1, we see that eq. (|5.35|) is still satisfied, while other 
matrix elements of E are transformed according to E — > RER T . One is free to choose R such that E in the new basis 
is diagonal. 

The above results pertain to a specific basis choice. With respect to an arbitrary basis, we showed in Section 3 that 
-Rcpi is a reflection through planes in K space, which implies that i?cpi is an improper rotation matrix that satisfies 

7? CP1 = 1 3 , Rcpi ¥= -la ■ (5-36) 

An explicit form is given by: 

R C pi = -R{n,ir), (5.37) 

where the unit vector h points in the direction normal to the reflection plane. In general, —R(h, it) possesses one 
non-degenerate eigenvalue —1 and a doubly-degenerate eigenvalue +1. We identify h as the eigenvector of —R(h, it) 
corresponding to the non-degenerate eigenvector —1. The eigenvectors corresponding to the doubly-degenerate eigen- 
values of i?cpi span the reflection plane in K space. 

The scalar potential exhibits the CP1 symmetry, eq. (|5.33|) . in some basis if 



K(x) ->• -R(h, vr)K(x') (5.38) 

is a symmetry. That is, the THDM parameters specified in eq. (|2.12l) must satisfy eq. (|3.40p with R = Rcpi, which 
yields: 



£ and if are eigenvectors of —R(n,ir) with eigenvalue +1 £ and rj are perpendicular to n, (5.39) 

ER(n, it) = R(n, n)E ^=^> Eh is parallel to h . (5.40) 

Following the derivation of eq. (|5.16j) . either En = or Eh is an eigenvector of —R(h,ir) with eigenvalue —1. Since 
the latter is non-degenerate, we again recover eq. (|5 . 16|) . It is a simple matter to confirm that eqs. (|5.16l) and (|5.39[) 
reduce to eq. (|5.35[) when h = y. 

We can introduce simultaneous eigenvectors of E and —R(h, it), denoted by rh and h X m, which satisfy eq. (|5.17|) 
[cf. footnote QT] . These two vectors correspond to the two-fold degenerate eigenvalue +1 of — _R(n, it). As in the 
treatment of the Z2 symmetry, the form of the CP1 symmetry transformation does not uniquely specify the basis 
in K space. To fix the basis completely, one must specify h and the eigenvectors of —R(h, tt) corresponding to the 
two- fold degenerate eigenvalue +1. 

In summary, the CP1 symmetry corresponds to K(ir) — > — i?(n, 7r)K(x') for some choice of h. Imposing this 
symmetry on the scalar potential requires that 

£ and 77 are perpendicular to h , (5-41) 
E is diagonal with respect to the basis {h , rh , h X rh} , (5-42) 
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where rh is a simultaneous eigenvector of —R(n, tt) and E, with — R(fi, 7r)m = rh. In this case, the choice of ft and 
rh uniquely fixes the basis in K space. 




• CP2 symmetry 

The CP2 symmetry class is generated by requiring that the scalar potential parameters specified in eq. (12.121) arc 
invariant under R = —R(n, 0) = — 1 3 corresponding to 28 = tt in eq. (|4.10|) . which yields 

R C P2 =%=\ 0-1 1. (5.43) 

V -v 

The form of i?cP2 is basis-independent, since R(U)Rcp2R T (U) — i?cP2 for any choice of U. We may now invoke 
Theorem 2 from Section 4. From eq. (j4.11j) . we then obtain 

# 3 Mi=-l3- (5.44) 

Invariance under i?cP2 = — 13 requires £ = and rj — but leaves E arbitrary. We can then perform a basis transfor- 
mation to diagonalize E, which leaves £ and rj unchanged. Invariance under all three individual CP1 transformations 
i?3, R 2 , and R\, taken together, implies that 

£ = r) = and E = diag(/zi, fx 2 , Ms) • (5-45) 

Clearly, R3, R2, and R\, when separately applied, yield a stronger constraint than — 13 alone. But the three CP1 
symmetries are physically equivalent to — 13, since the scalar potentials that result from the two procedures are 
related by the change of basis that diagonalizes E. As shown in [46|, this equivalence may no longer hold once the 
Yukawa sector is taken into account. 




• CP3 symmetry 

It is convenient to choose a basis in which the CP3 symmetry class is generated as follows. We require the 
transformation R' 2 e 1 [defined in eq. (I4.10j) ] to be a symmetry, for any choice of angle such that 29i ^ (mod tt). 
Theorem 2 can again be employed to express R 2 g 1 as a product of three CP1 transformations. Eqs. (|4.10p and (14. lip 
then yield: 

^cos26>i -sin 26»i 0\ 

= R 3 R 2ei R Sl - (5.46) 

Requiring CP3 to be a symmetry implies that the THDM parameters specified in eq. (|2.12p must satisfy eq. (13.40)) 
with R = R 2dl , which yields: 

£ = r\ = , E = diag(^i , ^1 , ^3) . (5.47) 

It is straightforward to check that imposing the three CP1 symmetries i?3, R 2 g 1 and Rg ± separately is equivalent to 
requiring invariance under R 2Sl . Similarly to the case of the U(l) Peccei-Quinn symmetry, the invariance of the scalar 
potential under R 26l , for any single particular value of 29± 7^ (mod tt), implies invariance under R 2g for all values 
of 6. 

The above results pertain to a specific basis choice. With respect to an arbitrary basis, we showed in Section 3 that 
-RcP3 is an improper rotation matrix that satisfies i?cp 3 7^ I3 (i-e., -Rcp3 is not a pure reflection in K space). The 
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most general form for the CP3 transformation is given by i?cP3 = —R(n, 9) for < 8 < ir. Thus the scalar potential, 
in a basis specified by n, exhibits the CP3 symmetry if it is invariant under the transformation 

K(i) ->• -R(n, 6)K(x') , forO<6»<7r. (5.48) 

Note that we have excluded the case of 8 = 0, which corresponds to the CP2 symmetry transformation, and the 
case of 8 = ir, which corresponds to the CP1 symmetry transformation, both of which have already been treated 
above. When 8 ^ (mod 7r), the improper rotation matrix —R(n, 8) possesses three non-degenerate eigenvalues: — 1, 
— e l9 and —e~ 10 , and h is the normalized eigenvector of —R(n, 8) with eigenvalue —1. As in the analysis of the U(l) 
Peccei-Quinn symmetry, it is convenient to introduce normalized eigenvectors of —R(fi,0), denoted by rh and rh* , 
corresponding to the eigenvalues — e l9 and —e~ l9 , respectively. 

Under the symmetry governed by eq. ()5.48[) . the THDM parameters specified in eq. (|2.12j) must satisfy eq. p.40|) 
with R = —R(n,9). Since —R(h,8) has no eigenvalue +1 for 8 ^ (mod 7r), it follows that: 

-R{h,6)i = i, —R(fi, 6)r> = rj , £ = T7 = 0, (5.49) 

ER{fi, 8) — R(h, 9)E Ev is parallel to v for v = h , rh , rh* . (5.50) 

The derivation of eq. (|5.50[) is identical to the one given in the analysis of the U(l) Peccei-Quinn symmetry. Thus, 
we again recover the results of eqs. (|5.30l) and (|5.3ip . 

EijHj cx rii and Eijirij cx m, . (5.51) 

For the same reasons given below eq. (|5.31l) . there is (at least) a two-fold degeneracy among the eigenvalues of E. Since 
the constraints imposed by eqs. (|5.49l) and (|5.51[) are independent of the angle 8, the invariance of the scalar potential 
under — R(h, 8) for any value of 8 ^ (mod ir) yields the CP3 symmetry, which in turn implies the invariance of the 
scalar potential under —R(h, 8) for all values of 8. As expected, when h = z and rh = -^=(x — iy), eqs. (|5.49[) and 
(|5.51[) lead to the previous results given in eq. (|5.47l) for any value of 9 ^ (mod 7r). 

By virtue of eq. (I5.51[) . the vectors comprising the orthonormal set {n,m,m*} are simultaneous eigenvectors of 
—R(h, 8) and E. This means that E is diagonal with respect to the basis {n , rh , rh*}. Since rh and rh* are uniquely 
determined by h [cf. eqs. (|A.13|) and (|A.14[) ] when 9 ^ (mod 7r), it follows that n uniquely fixes the basis in K 
space. 

• The SO(3)-symmetric scalar potential 

Suppose that the THDM scalar potential is invariant under all possible Higgs-family transformations [cf. eq. (|2.5[) ]. 
The U(2) invariance group contains a gauged hypercharge U(l) subgroup, which is automatically satisfied by the 
Higgs Lagrangian. The ungauged invariance group is therefore isomorphic to SO(3). The elements of this SO(3) 
group are in one-to-one correspondence with the matrices R(U) defined in eq. (|2.9I) . Hence, the invariance of the 
scalar potential under SO(3) requires that eq. (|2.13p must be satisfied for any SO(3) matrix R(U) £ SO(3). This 
clearly implies that 

£ = 0, r/ = 0, E = fil 3 , (5.52) 

which of course must hold for any choice of basis (since the set of all basis transformations has been promoted to a 
global symmetry of the scalar potential) . In K space, one rotation matrix is not sufficient to generate all the elements 
of the SO(3) symmetry class. Since any rotation matrix must be of the form R(h,9), the results of Section 5.1 
imply that a single rotation matrix will generate either the Z2 symmetry class (if 9 = n) or the U(l) Peccei-Quinn 
symmetry class, if < 9 < it. However, one can generate all the elements of the SO(3) symmetry class with two 
rotation matrices, R(hi 7 9) and R(h 2 , 9), where hi X n 2 ^ 0. In group theory language, this is just the statement that 
starting from two non-trivial rotations for which the corresponding rotation axes are non-collinear, it is possible to 
generate any element of SO(3) by successive multiplication of the two initial rotation matrices and products thereof. 
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For example, the S0(3) symmetry class can be generated through the application of both 

''cos «i 
R ai — sinai 

V 




with < ai < 7r, 



(5.53) 



and 



R'. 




With < «2 < 7T. 



(5.54) 



Since there is a infinite number of axes non-collinear with the axis of the first rotation, we have infinitely many ways 
to generate S0(3) through combinations of Higgs-family transformations. 

In fact, the maximal invariance group [orthogonal to the gauged hypercharge U(l)] of the S0(3)-symmetric scalar 



potential is 0(3), as noted in 4J, [50j. This can be seen either by examining the symmetry transformation properties 



of the scalar fields or via the transformations of the scalar field bilinears. Eq. (|5.52p implies that the scalar potential 
is invariant with respect to Higgs family transformations, tpi(x) — > Uijtpj(x), for any choice of the unitary matrix U. 
Since all potentially complex parameters of the S0(3)-invariant scalar potential are zero, eq. (|5.52p also implies 
that the scalar potential is invariant with respect to GCP transformations, <fi(x) — > Xijip*(x'), for any choice of 
the unitary matrix X. Mathematically, complex conjugation is isomorphic to the discrete group Z2, and hence the 
relevant invariance group is SO(3)<g>Z2 = 0(3). In K space, it immediately follows from eq. (|5 .52[) that eqs. (|2.13l) 
and (I3.40[) are satisfied for all proper and improper rotation matrices R and R, respectively. In particular, using the 
results of Section 5.1, we immediately conclude that the K space symmetries of the S0(3)-invariant scalar potential 
can also be generated by a product of two non-collinear CP3 transformations, —R(fii,d) and —R('h,2 ) 0), or by a 
product of the Higgs family transformation R(hi, 9) and a CP3 transformation —R(fi2, 0), where fii X 7^ 0. Thus, 
the full 0(3) group consisting of all proper and improper rotation matrices is generated. 14 



5.2. Identifying the THDM symmetry classes 

We have examined the constraints on the scalar potential parameters for each of the six symmetry classes listed 
in Table I. We can reverse the process and determine the symmetry class given the scalar potential parameters. We 
summarize the results of Section 5.1 in Table II, which agrees, of course, with the basis-independent characterization 



of the symmetry classes given in [391 . 144L 65]. In particular, given a set of constraints on £, 77 and the eigenvalues 
and eigenvectors of E, one can uniquely determine the maximal symmetry class of the scalar potential. Note that 
£ X 77 = for all symmetry classes, with the possible exception of CP1. In the case of CP1, the only necessary 
conditions on the scalar potential parameters are: 15 

£ • e, = 77 • e 4 = , (£,77)^(0,0), (5.55) 

where is one of the three eigenvectors of E (independently of the degeneracy of the eigenvalues of E). When the 
scalar potential parameters are generic, £ X 77 7^ 0. In this case the constraints given in eq. (|5.55[) are equivalent to 
the statement that £ X 77 is an eigenvector of E 39]. 16 On the other hand, if £ X 77 — 0, then the symmetry class 



The reader is cautioned that after the inclusion of the Yukawa sector, the corresponding symmetry properties of each of these pairs of 
transformations might be different from the pure scalar theory. 

Using the results of Table II, it is easy to check that £ • e = r\ ■ e = (where e is one of the eigenvectors of E) must hold for all symmetry 
classes. For the higher symmetry classes, additional constraints must be satisfied as exhibited in Table II. 

In this case, one can define n as the unit vector that points in the direction of £ X rj, and eq. Q5.41j l is automatically satisfied. 
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symmetry class constraints on £ and 77 constraints on E 

1*2 £ x e = 7] x e = 0; (£, rj) 7^ (0, 0) eigenvalues are non-degenerate 

or or 

£; X e — 77 X e — 0; (£, 77) 7^ (0, 0) E possesses doubly-degenerate eigenvalues 

U(l) £ X e — rj X e = 0; (£, T7) 7^ (0, 0) i? possesses doubly-degenerate eigenvalues 

or or 

£; X 77 = 0; (£ , 77) 7^ (0, 0) -B possesses triply-degenerate eigenvalues (E = /jj.3 

SO(3) {iiV) = (0,0) £ possesses triply-degenerate eigenvalues (E = /JI3 

CP1 £ X 77 is an eigenvector of E eigenvalues are unconstrained 

or or 
£ X 77 = 0; £ ■ e = 77 • e = 0; (£, tj) 7^ (0, 0); eigenvalues are non-degenerate 
neither £ nor 77 is an eigenvector of E 

or or 
£ X 77 = 0; £ • e — 77 ■ e' = 0; (£, 77) 7^ (0, 0); £ possesses doubly-degenerate eigenvalues 
neither £ nor 77 is an eigenvector of E 



CP2 


(€,»?) 


= (0,0) 


eigenvalues are non-degenerate 


CP3 


K.i?) 


= (0,0) 


i? possesses doubly-degenerate eigenvalues 



TABLE II: The symmetry classes and the corresponding constraints on the scalar potential parameters. The unit vector e is 
one of the three eigenvectors of E corresponding to a non-degenerate eigenvalue of E, and the unit vector e' is an eigenvector 
of E corresponding to a doubly-degenerate eigenvalue of E. Note that if (£, 77) = (0, 0), then the eigenvectors of E play no role 
in constraining the scalar potential parameters. For the U(l) symmetry class, e is uniquely identified as the eigenvector of E 
corresponding to the non-degenerate eigenvalue (if it exists). For the CP1 symmetry class, £ and 77 are necessarily orthogonal 
to one of the eigenvectors of E. If £ and 77 are collinear [and (£, 77) 7^ (0, 0)], then we must impose one additional condition to 
ensure that the constraints relevant for the Z2 or U(l) symmetry classes are not satisfied. 

of the scalar potential will be CP1 if and only if the constraints of the other symmetry classes listed in Table II are 
not satisfied. For example, consider the case where one eigenvalue of E is non-degenerate and corresponds to the 
eigenvector ej in eq. (|5.55[) . and the other eigenvalue of E is doubly-degenerate. Then eq. (j5.55[) necessarily implies 
that the collinear vectors £ and rj are also collinear with some linear combination of the other two eigenvectors ej 
U *) °f Ei an d the scalar potential therefore exhibits either the Z2 or U(l) symmetry. If the eigenvalues of E 
are non-degenerate, or if in eq. (|5.55|) corresponds to a doubly-degenerate eigenvalue of E, then it is possible that 
neither £ nor rj is an eigenvector of E, in which case the symmetry class of the scalar potential is CP1. 

Suppose that the scalar potential parameters satisfy the constraints of a particular symmetry class. By adding 
additional constraints, it is often possible to promote the scalar potential to a class with a larger symmetry. That is, 
there is a hierarchy of symmetries that can schematically be represented by the following chain: 

CP1 < Z 2 < I ^ I < CP3 < SO(3) . (5.56) 

The meaning of the above equation is the following. A scalar potential that exhibits a Z2 symmetry also must exhibit 
a CP1 symmetry. If the scalar potential exhibits a U(l) symmetry, then it must also exhibit the symmetries that 
precede it in the chain, and so on. Note that the U(l) symmetry does not imply a CP2 symmetry and vice versa. 
However, a CP3-symmctric potential must exhibit both the U(l) symmetry and the CP2 symmetry. Eq. (|5.56|) is 
a consequence of the results of Table [TTJ For example, if we take a CP2-symmetric scalar potential and impose the 
additional condition that one of the eigenvalues of E is doubly-degenerate, then we promote the symmetry to CP3. 
Likewise, if we take a U(l)-symmetric scalar potential (assuming E 7^ /1I3) and then impose the additional constraint 
£ = 77 = 0, then we promote the symmetry to CP3. The rest of the hierarchy of symmetries can be deduced in a 
similar manner. 

Note that there are cases in which the symmetry is not enhanced if the degeneracy of one of the eigenvalues of E 
is increased. Indeed, these cases do not correspond to new symemtry classes beyond the six listed in Table [TTJ This 
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conclusion is a consequence of a classification theorem, originally proved in 44| based on a geometrical analysis, 17 
which states that the only possible symmetry classes of the THDM are those listed in Table HH This result can be 
checked by employing the renormalization group equations (RGEs) for the scalar potential parameters. In particular, 
the constraints on the scalar potential parameters for a given symmetry class must be respected at any scale and 
hence invariant under renormalization group evolution. 

The RGEs for the scalar potential parameters in the general THDM can be found in |18j, 



48 



GG 



67| . One can check 



that all the parameter constraints exhibited in Table II are invariant under renormalization group running. However, 
if there are further accidental degeneracies among the eigenvalues of E, these degeneracies will not in general be 
renormalization group invariant. As a simple example, consider a U(l)-invariant scalar potential in a basis where 
A5 = A6 = A7 = 0. That is, we have from eq. f)2 . 12[) 



m? 2 



77oo = !(Ai+A 2 ) + iA 3 , 
E = i diag(A 4 , A 4 , |(Ai + A 2 ) - A 3 ) . 



\m = |(Ai- A 2 )y 



(5.57) 



In this case e = (0,0,1), and the constraints listed in Table II are clearly satisfied. Now, suppose we add the 
following additional constraint: Ai + A 2 = 2(A 3 + A4). In this case, the matrix E possesses one triply-degenerate 
eigenvalue, and E — /il 3 with fi = A4/4. However, this is not a new symmetry class because the constraint that 
Ai + A 2 = 2(A 3 + A4) is not in general renormalization group invariant. Using the RGEs for the scalar potential 
parameters with A5 — Ag = A7 = 0, we find: 



8tt 2 /3 Ai+A2 _ 2(A3+A4) = [Ai + A 2 - 2(A 3 + A 4 )] [3(A X + A 2 ) + 2A 4 ] + 3(A X - A 2 ) 2 



(5.58) 



Thus, if Ai =^ A 2 , then the condition Ai + A 2 — 2(A 3 + A4) — is not stable under renormalization group running. If 
in addition Ai = A 2 and m\ x — m 22 , then the scalar potential is SO(3)-invariant, corresponding to one of the other 
symmetry classes of Table II. Note that it is sufficient to set Ai = A 2 in eq. (|5.58[) to obtain a fixed point in the RGE 
at Ai + A 2 = 2(A 3 + A4) = 0. This would correspond to a softly broken SO(3)-invariant scalar potential, where the 
soft-breaking is due to m 2 ^ 7^ m\ 2 - Of course, this soft symmetry breaking is invisible to the RGE running of the 
coefficients of dimension-four terms of the Higgs Lagrangian. 

In general, additional degeneracies in the eigenvalues of E in Table II are either unstable with respect to renormal- 
ization group running or else correspond to a known symmetry class of enhanced symmetry that is either exact or 
softly-broken by dimension-two terms of the scalar potential. Thus, we conclude that no additional symmetry classes 
beyond those listed in Table II exist, in agreement with the classification theorem of |44j. 



5.3. Symmetries in the .E-diagonal basis 

The parameters of the scalar potential given by eq. (|2.2j) are described in the K space formalism by eqs. (12. lip 
and (I2.12[) . In particular, under U(2) basis transformations of eq. (|2.5p £ and 77 transform as real three- vectors 
and E transforms as a real Cartesian second-rank symmetric tensor under SO (3), which is related to the U(2) basis 
transformation of the scalar fields via eq. (|2.9[) . Starting in a generic basis, one can always transform to a special 
basis in which E is diagonal, 

RdERJ) = Ed = diag(^i, /x 2 , ^t 3 ) , (5.59) 



Note that in [44| the derivation of the symmetry classes assumes strong stability of the potential, that is, the quartic terms alone 
guarantee stability. Here we consider the general case including weak and marginal stability where the stability of the potential is 
guaranteed only after inclusion of the quadratic terms. In fact, the classification of Table ITT1 is valid even for unstable potentials. 
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rotation matrices in the i5-diagonal basis symmetry class 



constraints on scalar potential parameters 




$Xe — r/Xe — 0; eigenvalues of E are unconstrained 
£xn = 77Xn = 0; eigenvalues of E are degenerate 
£xe = 77Xe = 0; eigenvalues of E are doubly-degenerate 
£Xn = r]Xn = 0; E = fil 3 

£ ■ e — r) ■ e = 0; eigenvalues of E are unconstrained 
£ • n = r) ■ n — 0; eigenvalues of E are degenerate 
£ = T) = 0; eigenvalues of E are unconstrained 
£ = r\ = 0; eigenvalues of E are doubly-degenerate 
£ = ry = 0; E = ^1 3 



-R(n, 9), n e , < < m CP3 [SO(3)] 



TABLE III: The proper [and improper] rotation matrices R [and R] that generate the symmetry classes in K space, where the 
scalar fields are defined in the _B-diagonal basis. The unit vector e is any one of the three eigenvectors {ei, e.2, £3} of E. In the 
cases of the Z2 and CPl symmetry classes where ft/ej (for i = 1,2,3), E possesses doubly-degenerate eigenvalues if n can be 
expressed as a linear combination of two of the eigenvectors of E; otherwise E = /XI3. The last line of the table corresponds 
to a special case of the CP3 symmetry class in which all three eigenvalues of E are degenerate. The resulting scalar potential 
parameter constraints yield an SO(3)-invariant scalar potential. 

where the fa are the eigenvalues of E and Rd € SO(3). We begin by assuming that the scalar potential parameters 
are generic. In this case, the eigenvalues of E are non-degenerate, and Rd is unique up to the ordering of its rows 
and/or columns, which corresponds to a reordering of the diagonal elements of Ed- The corresponding normalized 
eigenvectors E are unique (up to an irrelevant multiplicative phase) and will be denoted e±, e.2 and e.3. Working in 
the basis where E is diagonal (henceforth called the i?-diagonal basis) is equivalent to expressing the matrix E with 
respect to the orthonormal basis {ei , e 2 , e 3 }. 

Having chosen the -E-diagonal basis, we now investigate the form of the various THDM symmetry classes. In each 
case, the symmetry transformation is of the form: 



where R is a proper rotation and R is an improper rotation. In Table III, we summarize the possible forms for 
R and R and identify the relevant symmetry class. We also indicate the corresponding constraints on the scalar 
potential parameters. In some cases, the imposition of the symmetry will require that the eigenvalues of E exhibit 
some degeneracy. If the symmetry imposes no such constraint, we say that the eigenvalues of E are unconstrained 
(i.e. non-degenerate for generic choices of the parameters of the scalar potential). 

We now discuss some of the salient points of Table III. First, the constraints on £ and r\ are precisely the same as the 
ones given in Table II. These constraints are determined simply by identifying the eigenvector of R or R (if it exists) 
corresponding to the eigenvalue +1. Indeed, when no eigenvalue +1 exists (as in the cases of CP2 and CP3), it follows 
that £ = 77 = 0. Second, in the cases of proper and improper rotation matrices parameterized in terms of n = e, where 
e is one of the eigenvectors of E, no constraints on the eigenvectors of E arise for the Z2, CPl and CP2 symmetry 
classes (in which case e can be any one of the three eigenvectors of E). For the U(l) and CP3 symmetry classes, the 
eigenvalues of R{h, 9) and R = —R(n, 9) are non-degenerate, and the corresponding eigenvectors orthogonal to n are 
complex conjugates of each other. Since E commutes with R, the eigenvectors of R are also eigenvectors of E. Finally, 
as the eigenvalues of the real symmetric matrix E must be real, it follows that the eigenvalue of E corresponding to 
the complex conjugate pair of eigenvectors is doubly-degenerate. 

However, if n 7^ (for i = 1, 2, 3), then additional constraints on the eigenvalues of E are imposed. Repeating the 
arguments of Section 5.1, we see that n must also be an eigenvector of E. This can only be consistent with h =/= ei 
if some of the eigenvalues of E are degenerate, in which case linear combinations of the within the degenerate 
subspace are also eigenvectors of E. For example, in the case of the Z2 or the CPl symmetry classes, if n is a 
linear combination of two of the eigenvectors of E, then eq. (|5 . yields a two-fold degeneracy in the eigenvalues 



K(x) -> iZK(ss) , 
K(ar) -> RK(x') , 



for Higgs family symmetries , 
for GCP symmetries , 



(5.60) 
(5.61) 
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of E. Likewise, if ft is a linear combination of all three of the originally chosen eigenvectors of E, then eq. (I5.16P 
yields a three-fold degeneracy in the eigenvalues of E, which means that E is proportional to the identity matrix 
(i.e. E = fits). Likewise, in the case of the U(l) or the CP3 symmetry classes, if n ^ e, (for i = 1,2,3), then the 
eigenvectors of R or R, denoted by rh and rh* in Section 5.1, are linear combinations of all three of the eigenvectors of 
E [cf. eq. (|A.13[) ]. Consequently, eq. (|5.31[) requires that all three eigenvalues of E must be degenerate, i.e. E = /Z.I3. 
Here is a simple example to illustrate the last point. Suppose we work in a coordinate system in which the eigenvectors 
of E are ei = (1, 0, 0), e.2 = (0, 1, 0), and £3 = (0, 0, 1). One possible choice for the axis of the rotation matrix R(h, 9) 
is h = ^7j(ei + e 2 ) = ^(1, 1,0). In this case, we use eq. (IA.13[) to obtain 

m = U^=(l + i) , (1 + *), + (5-62) 



2VV2 V2 

where rh is the eigenvector of R(n, 9) corresponding to the non-degenerate eigenvalue e l9 of R(n, 9) [under the 
assumption that < 9 < 71"]. When we impose the symmetry constraint, ER(h,9) = RE(n,9), it follows that rh 
must also be an eigenvector of E. This is only possible if E = as indicated in the last row of Table III. Any 
other choice of ft ^ e, (for i = 1,2, 3) would lead to the same conclusion. 

To summarize, we have examined possible symmetry transformations in the .E-diagonal basis, in which the vectors 
£ and rj and the matrix E are given with respect to the orthonormal basis constructed from the eigenvectors of 
E, namely {ei , (=2,63}. Applying a symmetry transformation ±R(h,9), where n ^ (for i — 1,2,3), yields the 
symmetry classes Z2, U(l), CP1, CP2 and CP3, but with extra degeneracies among the eigenvalues of E. Nevertheless, 
these extra degeneracies do not correspond in general to new enhanced symmetry classes as argued in Section 5.2. 
In particular, symmetry transformations defined with respect to a specific basis choice can yield constraints on the 
scalar potential parameters that are not stable with respect to renormalization group running. 

Of course, one can always redefine the scalar fields of the tree-level scalar potential to achieve a particular choice 
of basis. However, once higher order corrections are taken into account, the energy scale of the scalar potential 
parameters becomes relevant. For example, the diagonalization required to achieve the E-diagonal basis depends on 
the renormalization scale M. In general, the eigenvalues and the directions of the eigenvectors of E will not remain 
fixed under a change in the renormalization scale. Thus, parameter constraints defined in some specific basis relative 
to some scale M will not automatically be preserved at some other scale M' ^ M. When discussing restrictions 
on parameters of the tree-level scalar potential due to possible symmetries, one has to check explicitly that they are 
preserved by renormalization group evolution. 

As an example, we consider the GCP transformation 




for some fixed U) with < u> < tt/2, which corresponds to an improper rotation in K space given by 

/ cos 2uj sin 2ui 0^ 



(5.63) 



Rui = 



sin2w -cos2w | . (5.64) 
\ 1 



We require this to be a symmetry. Using eqs. (|A.8[) . (|A.10[) and (|A.11|) . it follows that 

R^ = ~R(n, it) , with n = (sincj , — cos to , 0) , (5.65) 

which means that R u generates a CP1 transformation. If we apply this transformation in the E-diagonal basis, where 
e.\ = (1, 0, 0), e.2 = (0, 1, 0), and £.3 — (0, 0, 1), then h = sinw e.\ — cosw e 2 must also be an eigenvector of E. Hence 
in the E-diagonal basis, E = diag(/ii, /13) possesses doubly-degenerate eigenvalues. 

However, the existence of doubly-degenerate eigenvalues is not stable under renormalization group evolution. To 
verify this statement, we shall work in a real basis (i.e. a basis in which all scalar potential parameters exhibited 
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in eq. (|2.2[) are real), which is stable under renormalization group running. Using eq. (|2.12l) . the matrix E then takes 
the form: 

A 6 - A 7 

A4-A5 I (.-..()(>) 

^(A!+A 2 )-A 3/ 

where A5, A6 and A7 are real. In the -E-diagonal basis, \§ — A7. Applying the symmetry constraints, eq. (|3.40[) . 
imposed by R u , yields A5 = and we see that E possesses a doubly-degenerate eigenvalue, A4. We first note that the 
diagonal form for E is not preserved under renormalization group running. In particular, in the real basis, 

8^ 2 /3 A6 -A 7 = (A 6 - A 7 ) [3(A x + A 2 ) + 2A 4 + 4A 5 ] + 3(A 6 + A 7 )(A X - A 2 ) . (5.67) 

Since Ai 7^ A2 in general, we see that the diagonal form for E is not stable. However, to check whether the CP1 
symmetry is enhanced, one must examine the eigenvalues of eq. (|5.66[) to see whether the doubly-degenerate eigenvalue 
is preserved or not under the renormalization group running. A straightforward computation shows that eq. (|5.66j) 
has a doubly-degenerate eigenvalue if 

D = A 5 (Ai + A 2 - 2A 3 - 2A 4 + 2A 5 ) - (A 6 - A 7 ) 2 = . (5.68) 

The beta function for D is given by 

Pd = (Ai + A 2 — 2A 3 — 2A 4 + 2\ 5 )j3\ 5 + A5/? Al+A2 _2A3-2A4+2A 5 — 2(A 6 — A 7 )/?a 6 -a 7 

= A 5 (/3 Al + /3a 2 - 2/?a 3 ~ 2/3 Al ) + (Ai + A 2 - 2A 3 - 2A 4 + 4A 5 )/3a 5 - 2(A 6 - A 7 )(/? Ae - P\ T ) . (5.69) 

Inserting the corresponding RGEs in the real basis 3, 48[ 66, 67 1 yields: 



8tt 2 Pd = 4D(A 1 +A 2 + A 3 + 2A4-A5) + 3(A 6 + A 7 ) 2 (Ai+A 2 -2A3-2A4 + 2A 5 ) + 3(A 1 -A 2 ) [A 5 (A x - A 2 ) - 2A 2 + 2A 2 ] . 

(5.70) 

Indeed, as long as Ai = A 2 and A 7 = — A6 is not satisfied, we see that D — is not a fixed point of eq. (|5.70[) . We 
recognize Ai = A 2 and A 7 = — A6 as the constraints on the dimensionless parameters of the scalar potential in the 
exceptional region of the parameter space (ERPS) identified in [5^]. When we impose these ERPS conditions, the 
CP1 symmetry is promoted to a CP3 symmetry 18 (modulo possible soft-symmetry breaking squared- mass terms). 
Outside of the ERPS, the double-degeneracy of the eigenvalue of E is not a renormalization group invariant, and the 
CP1 symmetry is not enhanced. 

In Table III, when h ^ e-i (for i = 1,2,3), the enhanced degeneracies in the eigenvalues of E are not stable with 
respect to renormalization group running in the cases of the Z 2 , U(l) and CP1 symmetry classes. However, in the case 
of the CP3 symmetry class, the constraints on the scalar potential parameters coincide with eq. (|5.52[) . which are the 
same constraints imposed by the SO(3)-invariant scalar potential. In this case, the resulting parameter constraints 
are invariant with respect to renormalization group running. That is, the CPS symmetry class with an enhanced 
degeneracy of the eigenvalues of E is promoted to the SO(3) symmetry class. In particular, it is possible to generate 
an SO(3) Higgs family symmetry by applying any particular CP3 transformation with n ^ ej (for i = 1,2, 3) in the 
-E-diagonal basis. This seems to be in conflict with the results of Section 5.1, where it was shown that two non-collinear 
rotation matrices are necessary in order to generate the SO (3) symmetry class in a generic basis. However, given that 
an independent rotation matrix Rd is required [cf. eq. (|5 . 59[) ] in order to define the -E-diagonal basis, it is perhaps 
not surprising that the SO(3) symmetry class can be generated by applying a single (improper) rotation matrix in 
the E-diagonal basis. 



In the .E-diagonal basis, the ERPS conditions on the dimensionless couplings of the scalar potential corresponds to Ai = A2 and 
^5 = ^6 = -V = 0- One is free to make a change of basis that simultaneously interchanges the rows and columns of E while keeping E 
diagonal. In the new basis, Ai = A2 and As = A7 = are maintained, while A5 = is transformed to A5 = Ai — A3 — A4. The latter 
reproduces the conditions for a CP3-symmetric scalar potential given in Table I of [5dl| . 
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6 CONCLUSIONS 



The classification of the possible THDM symmetries is best done in a generic basis, where the structure of the 
various symmetry classes is transparent, and the resulting constraints on the scalar potential parameters can be 
obtained that are covariant with respect to basis transformations. Although it is possible to perform an analysis of 
symmetry classes in a specific basis (the .E-diagonal basis is a convenient choice to consider for this purpose), the 
resulting classification is complicated by enhanced symmetry constraints that are typically not renormalization group 
invariant. Such enhanced symmetry points in the scalar potential parameter space are accidental in nature and are 
not indicative of any new symmetry structures. 



CONCLUSIONS 



It is known that there are only six classes of symmetry-constrained potentials in the THDM [44 ] . Specific im- 
plementations were later explored in [5(J, using Higgs-family and generalized CP (GCP) symmetries of the THDM 
potential. In this paper, we have presented an analysis of the symmetry classes, which applies to completely general 
scalar potentials. That is, our analysis applies to scalar potentials that are stable in the strong, weak, or marginal 
sense, or even unstable. Furthermore, we have pursued the geometric K space interpretation of Higgs-family and 
GCP transformations 38,0]; the former are proper SO(3) rotations and the latter are improper rotations of the 



field bilinears in K space. We have constructed the relevant rotations and shown explicitly their effects on the Higgs 
scalar potential. This combines some known results with new ones into a unified scheme and sets the framework for 
our analysis. 

The following results have thus been obtained. We have clarified the relation of the classifications of GCP transfor- 



mations of [45( and [50|. We have given a simple geometric proof relating SO(3) rotations to two reflections through 
planes and improper rotations to one or three reflections through planes in K space. We have shown that any Higgs- 
family transformation can be considered as a product of two CP1 transformations and any GCP transformation is 
either a CP1 transformation or a product of three CP1 transformations. Based on this result we have provided a 
geometric interpretation of the surprising result presented in (5o| that all Higgs-family and GCP symmetries in the 
THDM can be generated from suitable CP1 symmetries. 
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Appendix A: Properties of 3 X 3 proper rotation matrices 

In this Appendix, we review the properties of 3 x 3 proper rotation matrices. Most of this material is standard and 
be found in 68 



The most general three-dimensional proper rotation is represented by an S0(3) matrix, R(n 1 9), whose form is 
uniquely specified by an axis of rotation, n, and a rotation angle 9. Conventionally, a positive rotation angle cor- 
responds to a counterclockwise rotation. The direction of the axis is determined by the right hand rule. Simple 
geometrical considerations imply that: 

R(n,8 + 2nk) = R(fi,6), k = 0, ±1 ± 2 . . . , (A.l) 

[R(n, 9)]- 1 = R(n, -6) = R(-n, 9) . (A.2) 
Combining these two results, it follows that 

R(n,2ir-8)=R(-n,6), (A.3) 

which implies that any three-dimensional rotation can be described by a counterclockwise rotation by 9 about an 
arbitrary axis n, where < 9 < it. 19 However, for 9 — tt eq. (|A.3|) yields: 

R(n,Tr) = R(-n,w), (A.4) 

which means that for the special case of 9 = 7r, i?(n,7r) and i?(— n, tt) represent the same rotation. Finally, if 9 = 0, 
then R(n, 0) = I3 is the identity operator, independently of the direction of ft. 
An explicit form for a general three-dimensional proper rotation is given by: 

Rij(fi, 9) = Sij cos 9 + n,rij(l — cos 9) — ti^n^ sin 9 . (A. 5) 

This result simplifies for the case of 9 = tt, 

Rij(h, tt) = 2n. l n J - . (A. 6) 

It is noteworthy that R(h, 9) is a symmetric matrix if and only if 9 = (mod tt). 

Given the SO(3) matrix R(n,9), one can determine the corresponding angle of rotation 9 and axis of rotation n. 
By taking the trace of eq. (jA.5|) . we immediately obtain 

Tr R{n, 8) = 1 + 2 cos 9 = 1 + e iB + e~ ie . (A.7) 

It immediately follows that 

cos0 = |(Tr R- 1), (A.8) 
which determines 9 uniquely in the convention that < 9 < tt. The axis of rotation is given by: 

n = 1 (#32 - i? 23 , R13 - R31 , R21 - R12) , for 9 f mod tt . (A.9) 

\J\6 — Tr R){1 + Tr R) 

For 9 = (mod tt), R(h, 9) is symmetric and cannot be determined from eq. (|A.9I) . In this case, eq. (|A.8|) determines 
whether cos0 = +1 or cos0 = —1. If cos (9 = +1, then Rij = Sij and the axis n is undefined. If cos9 = —1, then 
eq. (|A.6j) determines the direction of n up to an overall sign. That is, 

h is undetermined if 9 = , 

n = [e 1 ^\{l + R 11 ),e 2 ^\{l + R 22 ),e^\{l + R^ , if0 = 7r, (A.10) 



In the convention adopted here, the overall sign of ft is meaningful for 6 =i mod ir. 
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where the individual signs £j = ±1 are determined up to an overall sign via 20 

eje.j = —== v = , for fixed i ^ j , Ru ^ -1 , ^ -1 . (A.ll) 

The ambiguity of the overall sign of n sign is not significant, since i?(n, 7r) and R(—h, tt) represent the same rotation 
as noted above [cf. eq. (|A.4[) ]. 

Since the axis of rotation h is invariant under all rotations R(n,6), it follows that h is an eigenvector of R(h,8) 
with corresponding eigenvalue +1, 

R(h,8)h = h. (A.12) 

Moreover, since i?(n, 8) is also a unitary matrix, we know that its eigenvalues are pure phases. Combining eq. (|A.7[) 
with det R(h, 6) = 1, it follows that the individual eigenvalues of R(n, 8) must be 1, e and e~ tB . For ^ (mod 7r), 
these three eigenvalues are non-degenerate. In this case, it is convenient to introduce normalized eigenvectors rh and 
rh* , corresponding to the eigenvalues e l6 and e~ %6 , respectively. Note that whereas n is a real vector corresponding 
to the axis of rotation, rh and its complex conjugate are complex vectors. The vectors n, rh and rh* are mutually 
orthonormal with respect to the inner product of a complex vector space (e.g., n • rh — rh ■ rh ~ and rh ■ rh* = 1). 
Explicitly, we have 

1 / in 2 (ni-in 2 ) . mi(ni-m 2 ) , . A „ , A1 ,, 

m = — n 3 H — , -«n 3 — , -(m - m 2 ) , for n ^ -z , (A. 13) 

V2 V 1 + "3 1 + ^3 / 

up to an overall phase factor that can be fixed by convention. This form is not very useful as 713 — > — 1. However, we 
can use eq. (|A.2[) to obtain (up to an overall phase) 

1 ( in 2 {n 1 +in 2 ) . mi(ni+m 2 ) . \ c ~ _l ~ t \ i a\ 

m = —= \ -n 3 , -in-i H , n± + vn 2 I , for n 2; . (A. 14) 

V2 V l~n 3 l-n 3 J 

Clearly, the eigenvectors n, rh and rh* are independent of the value of the rotation angle 9. In numerical work, it is 
convenient to use eq. (|A.13|) for n 3 > and eq. (|A.14j) for n 3 < 0. One can check that for n ^ ±£, eqs. (|A.13j) and 
(|A.14I) are identical up to an irrelevant multiplicative overall phase. 

In the case of 8 — tt, R(h,n) possesses one non-degenerate eigenvalue +1 and two degenerate eigenvalues — 1. 
The former is associated with the axis of rotation [cf. eq. (|A.12|) ]. The eigenvectors corresponding to the degenerate 
eigenvalues, denoted below by rh\ and rh 2 = n X rh 1 , can be chosen to be real and orthonormal. A convenient 
choice is: 21 

m x = [n 3 + — , — , -m , m 2 = — , n 3 + — , -n 2 . (A.15) 

\ 1 + n 3 1 + n 3 J \l + n 3 l + n 3 J 

However, any other orthonormal pair of vectors constructed from linear combinations of rhi and rh 2 would be equally 
suitable. Finally, in the case of 8 = 0, R{n,ir) = I3 possesses three degenerate eigenvalues +1. 
Finally, we prove an important result that is needed in the text. Let R £ SO(3) such that 

n=Rh. (A.16) 

Then 

R(h',8) = RR(h,8)R T . (A.17) 

To prove eq. (|A.17[) . we first note that 

Tr R(n,8) = Tr [RR(h, 8) R T ] = Tr R(h,8) , (A.18) 



20 If Ru = —1, where i is a fixed index, then rii = 0, in which case the corresponding n is not well-defined. 

21 In this case, any problem involving ns = — 1 can be avoided simply by employing eq. IIA.4II . 
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using the cyclicity of the trace and RR T = I3. It follows from eq. (IA.7I) that the angle of rotation of R(n,9) and 
RR(n, 9)R T must be the same. Next, we use eq. (|A.12[) to determine the axis of rotation of RR(h, 9)Ri T , 

RR(n, 9)R T (Rh) = RR(h, 9)n = Rn , (A.19) 

which implies that Rn is an eigenvector of RR(n,9)R T with eigenvalue +1. If 6 ^ 0, then the eigenvalue +1 is 
non-degenerate, in which case Rn is the axis of rotation of R(h',9), and eqs. (|A.16[) and (|A.17[) are confirmed. If 
9 = 0, then eq. (|A.17j) is trivially satisfied. 

Specializing to the case of n = z, it then follows that: 

R{n,8) = RR{z,9)R T , where h = Rz . (A.20) 



[1] K. Nakamura et al. [Particle Data Group], "Review of particle physics," J. Phys. G 37, 075021 (2010). 
[2] J.F. Gunion, H.E. Haber, G.L. Kane and S. Dawson, The Higgs Hunter's Guide (Westview Press, Boulder, CO, 2000). 
[3] T.D. Lee, "A theory of spontaneous T violation," Phys. Rev. D 8, 1226 (1973). 
[4] T.D. Lee, "CP nonconservation and spontaneous symmetry breaking," Phys. Rept. 9, 143 (1974). 
[5] R.D. Peccei and H.R. Quinn, "CP conservation in the presence of instantons," Phys. Rev. Lett. 38, 1440 (1977). 
[6] R.D. Peccei and H.R. Quinn, "Constraints imposed by CP conservation in the presence of instantons," 
Phys. Rev. D 16, 1791 (1977). 

[7] P. Fayet, "A gauge theory of weak and electromagnetic interactions with spontaneous parity breaking," 
Nucl. Phys. B 78, 14 (1974); 

[8] P. Fayet, "Supergauge invariant extension of the Higgs mechanism and a model for the electron and its neutrino," 
Nucl. Phys. B 90, 104 (1975). 

[9] K. Inoue, A. Kakuto, H. Komatsu and S. Takeshita, "Aspects of grand unified models with softly broken supersym- 
metry," Prog. Theor. Phys. 68, 927 (1982) [Erratum-ibid. 70, 330 (1983)]; "Renormalization of supersymmetry breaking 
parameters revisited," Prog. Theor. Phys. 71, 413 (1984). 
[10] R.A. Flores and M. Sher, "Higgs masses in the standard, multi-Higgs and supersymmetric Models," 
Annals Phys. 148, 95 (1983). 

[11] J.F. Gunion and H.E. Haber, "Higgs bosons in supersymmetric models. I," Nucl. Phys. B 272, 1 (1986) 

[Erratum-ibid. B 402, 567 (1993)]. 
[12] N.G. Deshpande and E. Ma, "Pattern Of Symmetry Breaking With Two Higgs Doublets," Phys. Rev. D 18, 2574 (1978). 
[13] H. Georgi, "A Model Of Soft CP Violation," Hadronic J. 1, 155 (1978). 

[14] H.E. Haber, G.L. Kane and T. Sterling, "The Fermion Mass Scale And Possible Effects Of Higgs Bosons On Experimental 

Observables," Nucl. Phys. B 161, 493 (1979). 
[15] J.F. Donoghue and L. F. Li, "Properties Of Charged Higgs Bosons," Phys. Rev. D 19, 945 (1979). 

[16] E. Golowich and T.C. Yang, "Charged Higgs Bosons And Decays Of Heavy Flavored Mesons," Phys. Lett. B 80, 245 
(1979). 

[17] L.J. Hall and M.B. Wise, "Flavor Changing Higgs Boson Couplings," Nucl. Phys. B 187, 397 (1981). 

[18] H.E. Haber and R. Hempfling, "The Renormalization group improved Higgs sector of the minimal supersymmetric model," 

Phys. Rev. D 48, 4280 (1993) arXiv:hep-ph/9307201 . We alert the reader to one typographical error in Appendix A of 

this paper. In eq. (A2), replace 12Ag with I2A7. Note that the renormalization equations given in Appendix A were derived 

for the CP-conserving THDM in which all scalar potential parameters are real. 
[19] G. Cvetic, "CP violation in bosonic sector of SM with two Higgs doublets," Phys. Rev. D 48, 5280 (1993) 

[arXiv:hep ph/9309202] . 

[20] F.J. Botella and J. P. Silva, "Jarlskog-like invariants for theories with scalars and fermions," Phys. Rev. D 51, (1995) 3870 
[arXiv:hep-ph/9411288 . 

[21] L. Lavoura and J. P. Silva, "Fundamental CP violating quantities in a SU(2) x U(l) model with many Higgs doublets," 

Phys. Rev. D 50, 4619 (1994) [arXiv:hep-ph/9404276j . 
[22] L. Lavoura, "Signatures of discrete symmetries in the scalar sector," Phys. Rev. D 50, 7089 (1994) [arXiv:he p-ph/9405307 . 
[23] W. Bernreuther and O. Nachtmann, "Flavour dynamics with general scalar fields," Eur. Phys. J. C 9, 319 (1999) 

[arXiv:hep-ph/9812259] . 



29 



A PROPERTIES OF 3 x 3 PROPER ROTATION MATRICES 

I.F. Ginzburg and M. Krawczyk, "Symmetries of two Higgs doublet model and CP violation," 
Phys. Rev. D 72, 115013 (2005) [arXiv:hep-ph/0408Qlll . 

S. Davidson and H.E. Haber, "Basis-independent methods for the two-Higgs-doublet model," 
Phys. Rev. D 72, 035004 (2005) [Erratum-ibid. D 72, 099902 (2005)] [arXiv:hep-ph/0504050 . 



J.F. Gunion and H. E. Haber, "Conditions for CP-violation in the general two-Higgs-doublet model," 
Phys. Rev. D 72, 095002 (2005) [arXiv:hep-ph/0506227| . 

H.E. Haber and D. O'Neil, "Basis-independent methods for the two-Higgs-doublet model. II: The significance of tan 
Phys. Rev. D 74, 015018 (2006) [arXiv:hep-ph/0602242| . 

R. Barbieri and L.J. Hall, "Improved naturalness and the two Higgs doublet model," arXiv:hep-ph/0510243. 

G.C. Branco, M.N. Rebelo and J.I. Silva-Marcos, "CP-odd invariants in models with several Higgs doublets," 

Phys. Lett. B 614, 187 (2005) [arXiv:hep-ph/0502118l . 

A. Barroso, P.M. Ferreira and R. Santos, "Charge and CP symmetry breaking in two Higgs doublet models," 



Phys. Lett. B 632, 684 (2006) arXiv:hep-ph/0507224 



R. Barbieri, L.J. Hall and V.S. Rychkov, "Improved naturalness with a heavy Higgs: An alternative road to LHC physics,' 



Phys. Rev. D 74, 015007 (2006) arXiv:hep-ph/0603188 



L. Fromme, S.J. Huber and M. Seniuch, "Baryogenesis in the two-Higgs doublet model," JHEP 0611, 038 (2006) 
[arXiv:hep-ph/0605242j . 

A. Barroso, P.M. Ferreira and R. Santos, "Neutral minima in two-Higgs doublet models," Phys. Lett. B 652, 181 (2007) 
[arXiv:hep-ph/0702098] . 

J.M. Gerard and M. Herquet, "A twisted custodial symmetry in the two-Higgs-doublet model," 
Phys. Rev. Lett. 98, 251802 (2007) [arXiv:hep-ph/0703051| . 

F. Mahmoudi and O. Stal, "Flavor constraints on the two-Higgs-doublet model with general Yukawa couplings," 
[arXiv:0907. 17911 [hep-ph]]. 

P.M. Ferreira and J. P. Silva, "A Two-Higgs Doublet Model With Remarkable CP Properties", accepted for publication in 
Eur. Phys. J. C, [arXiv:1001.0"574V 2 [hep-ph]]. 

J. Velhinho, R. Santos and A. Barroso, "Tree level vacuum stability in two-Higss doublet models," 
Phys. Lett. B 322, 213 (1994). 

F. Nagel, "New aspects of gauge-boson couplings and the Higgs sector," Ph. D. thesis (University of Heidelberg, 2004), 
available from http://archiv.ub.uni-heidelberg.de/archiv/4803. 

LP. Ivanov, "Two-Higgs-doublet model from the group-theoretic perspective," Phys. Lett. B632, 360 (2006) 
[hep-ph/0507132l . 

M. Maniatis, A. von Manteuffel, O. Nachtmann and F. Nagel, "Stability and symmetry breaking in the general two-Higgs- 



doublet model," Eur. Phys. J. C 48, 805 (2006) arXiv:hep-ph/0605184 



M. Maniatis, A. von Manteuffel and O. Nachtmann, "Determining the global minimum of Higgs potentials via Groebner 
bases-applied to the NMSSM," Eur. Phys. J. C 49, 1067 (2007) [arXiv:hep-ph/0608314| . 

C.C. Nishi, "CP violation conditions in 7V-Higgs-doublet potentials," Phys. Rev. D 74, 036003 (2006) 



[Erratum-ibid. D 76, 119901 (2007)] arXiv:hep-ph/0605153 



LP. Ivanov, "Minkowski space structure of the Higgs potential in 2HDM," Phys. Rev. D 75, 035001 (2007) 
[Erratum-ibid. D 76, 039902 (2007)] [arXiv:hep-ph/0609018| . 

LP. Ivanov, "Minkowski space structure of the Higgs potential in 2HDM. II. Minima, symmetries, and topology," 
Phys. Rev. D77, 015017 (2008) |arXiv: 0710.34901 [hep-ph]]. 

M. Maniatis, A. von Manteuffel and O. Nachtmann, "CP Violation in the general two-Higgs-doublet Model: a geometric 
view," Eur. Phys. J. C 57, 719 (2008) [arXiv:0707.3344l [hep-ph]]. 

M. Maniatis, A. von Manteuffel and O. Nachtmann, "A new type of CP symmetry, family replication and fermion mass 
hierarchies," Eur. Phys. J. C 57, 739 (2008) [arXiv:0711.3760l [hep-ph]]. 

M. Maniatis and O. Nachtmann, "On the phenomenology of a two-Higgs-doublet model with maximal CP symmetry at 
the LHC," JHEP 0905, 028 (2009) [arXiv:0901.434"T1 [hep-ph]] . 

E. Ma and M. Maniatis, "Symbiotic Symmetries of the Two-Higgs-Doublet Model," Phys. Lett. B 683, 33 (2010) 
[arXiv:0909. 28551 [hep-ph]]. 

M. Maniatis and O. Nachtmann, "On the phenomenology of a two-Higgs-doublet model with maximal CP symmetry at 
the LHC, II: radiative effects," [arXiv:0912.2727l [hep-ph]]. 

P.M. Ferreira, H.E. Haber and J. P. Silva, "Generalized CP symmetries and special regions of parameter space in the 
two-Higgs-doublet model," Phys. Rev. D 79, 116004 (2009) [arXiv:0902. 15371 [hep-ph]] . 



30 



A PROPERTIES OF 3x3 PROPER ROTATION MATRICES 



[51] G. Ecker, W. Grimus and W. Konetschny, "Quark mass matrices in left-right symmetric gauge theories," 
Nucl. Phys. B 191, 465 (1981). 

[52] G. Ecker, W. Grimus and H. Neufeld, "Spontaneous CP violation in left-right symmetric gauge theories," 
Nucl. Phys. B 247, 70 (1984). 

[53] G. Ecker, W. Grimus and H. Neufeld, "A standard form for generalized CP transformations," J. Phys. A 20, L807 (1987). 
[54] H. Neufeld, W. Grimus and G. Ecker, "Generalized CP invariance, neutral flavor conservation and the structure of the 

mixing matrix," Int. J. Mod. Phys. A 3, 603 (1988). 
[55] G.C. Branco, L. Lavoura and J. P. Silva, CP Violation (Oxford University Press, Oxford, UK, 1999). 
[56] O. Nachtmann, Elementary Particle Physics: Concepts And Phenomena (Springer- Verlag, Berlin, Germany, 1990). 
[57] I.I. Bigi and A.I. Sanda, CP Violation, 2nd edition (Cambridge University Press, Cambridge, UK, 2009). 
[58] J. Bernabeu, G.C. Branco and M. Gronau, "CP restrictions on quark mass matrices," Phys. Lett. B 169, 243 (1986). 
[59] H.K. Dreiner, H.E. Haber and S.P. Martin, "Two-component spinor techniques and Feynman rules for quantum field theory 

and supersymmetry," Phys. Rept. 494, 1 (2010) [arXiv:0812.T594l [hep-ph]]. 
[60] R.A. Horn and V.V. Sergeichuk, "Canonical forms for unitary congruence and 'congruence," Linear and Multilinear 

Algebra 57, 777 (2009) [arXiv:0710.1530l [math.RT]] . 
[61] W. Grimus, G. Ecker, "Basis transformations in generation space and a criterion for the existence of standard forms for 

unitarily congruent matrices," J. Phys. A 21, 2825 (1988). 
[62] S.L. Glashow and S. Weinberg, "Natural conservation laws for neutral currents," Phys. Rev. D 15, 1958 (1977). 
[63] E.A. Paschos, "Diagonal neutral currents," Phys. Rev. D 15, 1966 (1977). 

[64] P.M. Ferreira and J. P. Silva, "Discrete and continuous symmetries in multi-Higgs-doublet models," 

Phys. Rev. D 78, 116007 (2008) [arXiv: 0809. 27881 [hep-ph]] . 
[65] P.M. Ferreira, M. Maniatis, O. Nachtmann and J. P. Silva, "CP properties of symmetry-constrained two-Higgs-doublet 

models," JHEP 1008, 125 (2010) [arXiv: 1004.32071 [hep-ph]]. 
[66] W. Grimus and L. Lavoura, "Renormalization of the neutrino mass operators in the multi-Higgs-doublet standard model," 

Eur. Phys. J. C 39, (2005) 219 [arXiv:hep-ph/04092311 . 
[67] P.M. Ferreira and D.R.T. Jones, "Bounds on scalar masses in two Higgs doublet models," JHEP 0908 (2009) 069 

|arXiv:0903.2"856l [hep-ph]] . 
[68] W.-K. Tung, Group Theory in Physics (World Scientific, Singapore, 1985). 

[69] K.N. Srinivasa Rao, The Rotation and Lorentz Groups and their Representations for Physicists (Wiley Eastern Limited, 
New Delhi, 1988). 

[70] S.L. Altmann, Rotations, Quaternions, and Double Groups (Dover Publications, Inc., New York, NY, 2005). 



31 



